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I .  INTRODUCTION 


The  majority  of  information  on  the  burning  behavior  of  propellants  used 
in  interior  ballistic  gun  simulations  is  obtained  through  burning  samples  of 
the  propellant  in  a  closed  bomb.  Analysis  of  the  resulting  pressure-time 
data,  either  in  terms  of  burning  rate  information  or  experimental  impetus,  is 
of  vital  importance  to  the  interior  ballistician  and  propellant  formulator. 
Highly  sophisticated  computer  programs  have  been  written  to  perform  the  needed 
burning  rate  analyses .  Results  from  these  computations  have  proven  to  be 
extremely  valuable  in  predicting  propellant  performance  in  guns.  Such  burning 
rate  data,  for  instance,  when  used  in  two-phase-flow  interior  ballistic  codes 
such  as  NOVA  have  permitted  prediction  of  not  only  gun  muzzle  velocities  and 
peak  chamber  pressures  but  pressure  oscillations  as  well.^ 

Unfortunately,  sophistication  is  obtained  at  a  price.  Often  the  length 
or  other  requirements  of  the  burning  rate  reduction  program  limit  the 
equipment  on  which  it  can  be  used.  In  addition,  program  complexity  often 
hinders  quick  turn  around  between  closed  bomb  firings  and  data  reductions . 
Finally,  the  intricate  nature  of  the  program  can  inhibit  attempts  at 
modification  or  extension.  This  final  point  is  especially  critical  in 
research  programs  where  some  parameters,  such  as  grain  geometry,  may  fall 
outside  the  "normal  range"  accounted  for  in  the  program. 

The  drawbacks  mentioned  above  have  become  apparent  during  the  past 
several  years  in  our  work  on  new  propellants  for  advanced  ballistic 
applications.  The  program  under  discussion  was  written  as  an  answer  to  these 
needs.  In  addition  to  burning  rate  computations,  the  program  also  calculates 
the  propellant  experimental  impetus  based  on  maximum  experimental  pressure,  a 
highly  useful  feature  in  evaluating  experimental  performance  of  new 
formulations  relative  to  their  predicted  thermodynamic  performance  via  BLAKE.5 

MINICB  is  offered  as  a  practical  tool  for  propellant  formulators  and 
combustion  researchers  to  furnish  quick  feedback  on  results  of  closed  bomb 


1 


firings  which  could  be  used  in  guiding  formulation  efforts  or  combustion 
diagnostic  work.  Specifically,  MINICB  was  developed  to: 


a.  be  simple  and  easy  to  operate. 

b.  be  machine  independent.  In  fact,  MINICB  will  run  on  any 
computer  using  BASIC  and  on  some  programmable  hand 
calculators . 

c.  be  easily  modified,  especially  in  terms  of  non-standard 
grain  geometry. 

d.  be  usable  with  a  full  set  of  pressure-time  data  or  a 
limited  number  of  pressure- time  points 

e.  be  sufficiently  accurate  so  as  to  provide  useful 
information. 

In  this  report  we  present  the  formulation  of  the  underlying  equations 
together  with  the  simplifying  assumptions  which  were  used.  Validation  of  the 
program  and  suggestions  for  implementation  are  also  included. 


II.  DERIVATION  OF  EQUATIONS 


BURNING  RATE  EQUATIONS 

A.  MASS  FRACTION  BURNT 


In  determining  burning  rates  a  sample  of  propellant  together  with  an 
appropriate  amount  of  r'gnieer  is  ignited  in  a  closed  chamber..  As  the  igniter 
and  propellant  burn  the  pressure  within  the  chamber  is  recorded,  using  the 
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pressure- time  history  from  the  closed  chamber,  physical  properties  of  the 
closed  chamber,  thermochemical  properties  of  the  igniter  and  propellant,  and 
the  grain  geometry  of  the  propellant,  a  relation  between  the  burning  rate  of 
the  propellant  and  pressure  can  be  determined. 

The  burning  rate  analysis  is  a  dynamic  process  in  which  calculations  are 
performed  stepwise  along  the  pressure- time  history  of  the  burn.  Mass  and 
energy  balance  must  be  accounted  for  in  the  conversion  of  solid  propellant  to 
gaseous  products.  The  process  is  not  adiabatic  since  heat  loss  occurs  to  the 
walls  of  the  chamber.  The  treatment  of  heat  loss  may  be  simplified,  however, 
by  assuming  that  it  is  a  constant  percentage  throughout  the  burning  cycle  and 
that  it  can  be  accounted  for  by  adjusting  the  theoretical  gas  temperature  of 
the  system  by  the  ratio  of  the  observed  to  theoretical  maximum  pressures. 

Summarizing,  the  assumptions  made  in  the  burning  rate  analysis  are 
presented  below. 


1.  The  original  volume  of  the  closed  chamber  is  known. 

2.  The  igniter  is  completely  burnt  before  the  propellant 
begins  to  burn. 

3.  The  maximum  pressure  from  the  firing  is  known. 

4.  All  thermochemical  properties  of  the  igniter  and 
propellant  are  known. 

5.  The  Noble-Able  equation  of  state  can  be  used  to  describe 
the  state  of  the  closed  chamber. 

6.  The  gas  products  of  the  igniter  and  propellant  are 
thoroughly  mixed  at  all  times . 

7.  The  rate  of  heat  loss  is  constant  throughout  the  burning 
process  and  can  be  represented  as  stated  above. 

8.  The  temperature  of  the  system  before  heat  loss  can  be 
given  as  a  weighted  average  of  the  theoretical  flame 
temperatures  of  the  igniter  and  propellant  gases. 
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The  burning  rate  analysis  requires  equations  to  describe  the  mass  of 
propellant  burned  in  the  chamber  at  any  given  instant  in  terms  of  measured 
pressure  and  known  parameters.  To  derive  the  required  equations  start  with 
the  Noble-Able  equation  of  state. 

P(Vb  -  b)  «-  NRT  (1) 

Since  the  process  is  dynamic,  the  free  volume  of  the  bomb,  Vb,  will  be 
constantly  changing  and  is  given  by  equation  (2) . 

Vb  “  V  -  (C  -  Cp)/pp  -  (Cis  -  C±)/?i  (2) 

Also,  the  co -volume  correction  term  is  given  by, 

b  -  n^  +  npCp  (3) 

Now,  recalling  the  assumption  that  the  igniter  is  all  burnt  at  the  beginning 
of  the  analysis,  hence,  -  Cis,  and  substituting  (2)  and  (3)  into  (1)  gives 
equation  (4) . 

P(V  -  c/pp  +  Cp/pp  -  Cptip  -  CiH)  -  NRT  (4) 

Assuming  a  well  mixed  combination  of  the  gases  from  the  igniter  and 
propellant, 


N  =  Nt  +  Np  (5) 

equation  (4)  becomes: 

P(V  -  C/Pp  +  Cp/Pp  -  Cptip  -  C^)  -  RT(Nt  +  Np)  (6) 

Next  using  the  relation,  number  of  moles  =  weight/molecular  weight,  to  modify 
equation  (6) . 
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P(V  -  c/pp  +  cp/pp  -  Cptip  -  C^)  -  RT(Cp/mp  +  Ci/mi) 


(7) 


Now  equation  (7)  can  be  solved  for  C 


P* 


P[V  -  C/Pp  -Cini]  -  (RT/mi)Ci 


(RT/mp)  +  P[np  -  l/pp] 


(8) 


The  term  T  in  equation  (8)  stands  for  the  system  gas  temperature.  In  an 
adiabatic  system  this  would  be  a  weighted  average  for  the  igniter  and 
propellant  flame  temperatures  with  suitable  provisions  for  differences  in 
their  respective  heat  capacities.  Assuming  equivalent  heat  capacities  for 
propellant  and  igniter  combustion  products,  it  can  be  approximated  as 


T 


CTp  +  CiTi 


C  +  Ci 


(9) 


Assuming  further  that  the  experimental  shortfall  of  pressure  is  entirely  due 
to  cooling  of  the  gaseous  combustion  products,  the  system  temperature  may  be 
approximated  as 


CTp  +  Pmax 

X  -  -  *  -  (10) 

C  +  Ptheo 
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For  systems  where  the  amount  of  propellant  is  much  greater  than  the  amount  of 
igniter,  T  may  be  approximated  throughout  as  a  constant  defined  by  the  all- 
burnt  condition  of  equation  (10).  Ptheo  is  computed  using  the  Noble-Able 
equation  of  state  with  the  igniter  and  propellant  all  burnt.  This  is  equation 
(11)  given  below. 


Ptheo 


RT^/n^  +  RTpC/mp 


(V  -  niGi  -  ripC) 


(11) 


One  last  simplification  can  be  made  to  several  of  the  equations  by  using  the 
following  expression  for  impetus  of  the  igniter. 

F.  =  RT i/mi  (12) 

Equation  (11)  can  be  rewritten  as 


Ptheo 


FiCi  +  RTpC/“p 


(V  -  n^  -  UpC) 


(13) 


and  equation  (8)  for  the  mass  of  propellant  burned  at  any  instant  becomes 


P[V  -  C/Pp  -C^]  -  (FiT/Ti)Ci 


(RT/mp)  +  P[np  -  l/pp] 


(14) 
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Having  derived  the  equation  for  the  mass  of  propellant  burned  at  any 
instant  in  terms  of  measured  pressure  and  known  variables,  several  approaches 
to  computing  burning  rate  are  possible.  These  are  all  predicated  on  our 
ability  to  further  describe  the  mass  of  propellant  burned  in  terms  of 
geometric  changes  in  the  burning  propellant  grains.  For  example,  if  x 
represents  the  distance  burned  and  A  the  surface  area  of  the  propellant  when  x 
is  the  distance  burned  then 

Cp  -  A  *  x  *  pp  (15) 

or 

x  -  Cp/A*pp  (16) 

Finally,  differentiating  both  sides  of  (19)  with  respect  to  time  gives  an 
expression  for  the  burning  rate,  r. 

r  -  dx/dt  -  (dCp/dt)/A*pp  (17) 

This  is  the  approach  which  will  be  followed  here  since  knowledge  of  the 
instantaneous  surface  area  is  a  valuable  tool  in  analyzing  and  understanding 
closed  bomb  results.  Alternately,  x,  the  distance  burned,  hence  burning  rate, 
could  be  expressed  in  terms  of  the  mass  fraction  of  propellant  burned  without 
invoking  surface  area  at  all.  That  technique,  particularly  suited  to  hand 
calculators,  is  covered  in  Appendix  E. 

B.  SURFACE  AREA 

In  computing  the  burning  rate,  r,  in  equation  (17)  the  value  for  the 
surface  area,  A,  must  be  known  at  each  time  step  the  calculations  are  being 
performed.  For  this  program  the  surface  area  is  calculated  using  the 
following  scheme . 
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1.  The  amount  of  propellant  burnt,  Cp,  is  equal  to  the 
density  times  the  volume  of  propellant  burnt. 

2.  The  volume  of  propellant  burnt  is  determined  from  the 
initial  dimensions  of  the  propellant  grain  and  the 
linear  depth  burnt,  which  is  denoted  by  x. 

3.  From  1  and  2,  Cp  is  a  function  of  initial  propellant 
grain  dimensions,  density,  and  depth  burnt. 

4.  Thus,  the  depth  burnt  can  be  solved  for  in  terms  of 
initial  grain  geometry,  amount  burnt,  and  density. 

5.  The  surface  area  at  any  time  is  a  function  of  initial 
grain  geometry  and  depth  burnt  at  that  time. 

To  illustrate  the  scheme,  the  formula  used  in  MINICB  for  the  surface  area 
of  a  single  perf  grain  will  be  derived.  The  derivation  of  other  surface  area 
equations  can  be  found  in  Appendix  A  and  in  the  text  by  Hunt .  ^  For  the 
derivation  let  the  propellant  density  be  given  by  p  instead  of  Pp. 

Figure  1.  shows  a  single  perforated  grain  after  it  has  burnt  a  linear 
distance,  x,  on  all  surfaces.  The  volume  used  or  burnt  is  given  in  equation 
(18). 


Volume  Used  -  L*Pi*[(D/2)2  -  ((D  -  2x)/2)2] 

+  L*Pi*[((d  +  2x)/2)2  -  (d/2)2] 

+  2x*Pi*[ ( (D  -  2x)/2)2  -  ( (d  +  2x)/2)2]  (18) 


In  equation  (18),  the  first  term  represents  the  volume  used  along  the  outer 
lateral  surface  and  the  second  term,  the  volume  along  the  lateral  perf 
surface.  The  final  term  is  the  additional  volume  lost  from  the  ends  not 
accounted  for  in  the  other  two  terms.  Another  approach  to  determining  the 
volume  used  would  be  to  take  the  original  volume  and  subtract  the  volume 
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Now  the  total  amount  of  propellant  burnt  can  be  expressed  in  the  following 
relation. 


Cp  =  p  *  Volume  Used  (20) 

Substituting  (19)  into  (20) 

Cp  =  p  *  L*Pi*[xD  -  x2]  +  L*Pi*[xd  +  x2] 

+(x*Pi/2)*[D2  -  d2  -  4xD  -4xd]  (21) 

At  this  point  equation  (21)  can  be  solved  for  x  which  then  can  be  used  to 
determine  the  surface  area.  Expanding  and  rearranging 

2*Cp/(L*p)  -  (-4D  -  4d)*x2  +  (-d2  +  D2  +  2LD  +  2Ld)*x  (22) 

This  is  now  a  quadratic  equation  in  x  and  can  be  solved  using  the  quadratic 
formula.  First,  rearrange  equation  (22). 

(4D  +  4d)*x2  -  (-d2  +  D2  +2LD  +2Ld)*x  +  2*Cp/(L*p)  -  0  (23) 

Let  AA  «■  4D  +  4d,  BB  =  (-d2  +  D2  +2LD  +2Ld) ,  and  CC  »  2*Cp/(L*p) ,  then 


X 


BB  +-  V(-BB)2  *  4*AA*CC 


(24) 


2*AA 


At  this  point  it  is  necessary  to  determine  which  root  corresponds  to  the 
correct  value  of  x.  Several  observations: 

] .  AA  and  CC  are  both  positive. 

2.  BB  is  positive  since  D  is  large  than  d.. 

3.,  The  radical  must  be  real  since  a  solution  must  exist. 

4.  The  radical  is  always  less  than  BB  since  4*AA*CC  is 
positive  and  is  subtracted  from  BB  . 

5.  Based  upon  1-4,  there  are  two  positive  real  roots  for 
x. 

6.  Now,  x  should  be  the  smaller  of  the  roots.  If  x  was  the 
larger  root  then  the  smaller  root  would  also  be  a  valid 
solution.  However,  this  would  result  in  two  possible 
solutions  which  is  physically  not  possible. 

Thus,  x  must  correspond  to  the  root  with  the  negative 
sign. 

Therefore , 


x 


4*AA*CC 


2*AA 


(25) 


Referring  to  Figure  1,  the  surface  area  of  the  grain  corresponding  to  a  depth 
x  burnt  is  given  in  equation  (26). 
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A  -  (L  -  2x)*Pi*[D  -  2x]  +  (L  -  2x)*Pi*[d  +  2x] 


+  2*Pi*[((D  -  2x)/2)2  -  ((d  +  2x)/2)2]  (26) 

This  equation  can  now  be  used  to  determine  the  surface  area  of  the  grain  for  a 
given  depth  burnt. 

Summarizing,  the  burning  rate  is  given  by  equation  (17) 

r  *=  dx/dt  -  (dC  /dt)/A*p  (17) 

r  r 

To  determine  the  rate  the  following  procedure  is  used. 

1.  From  equation  (14)  Cp  can  be  determined  for  each  value 
of  pressure  in  the  pressure-time  history. 

2 .  The  values  for  dCp/dt  can  be  determine  from  Cp  using 
a  numerical  technique. 

3.  Using  equation  (26)  or  a  similar  one  for  the  appropriate 
grain  geometry  the  values  for  the  instantaneous 
surface  area  can  be  determined  for  each  value  of  Cp 
determined  in  1. 

4.  Since  pp  is  a  constant,  the  burning  rate  can  now  be 
readily  calculate  from  equation  (17). 

EXPERIMENTAL  IMPETUS 

The  traditional  approach  for  determining  experimental  impetus  is 
presented  by  Hunt^  in  his  text  on  interior  ballistics .  Exact  details  can  be 
found  in  Appendix  D.  This  method  requires  a  series  of  closed  bomb  firings 
with  different  loading  densities.  In  some  propellant  research  efforts,  where 
burning  rate  is  of  primary  interest  and  sample  amounts  are  limited,  it  is 


often  not  practical  to  fire  such  a  series.  If  a  reasonably  accurate  value  of 
the  propellant  co-volume  can  be  obtained,  however,  then  a  reliable  measure  of 
experimental  impetus  may  be  obtained  with  a  single  closed  bomb  firing. 

In  propellant  development  efforts  it  is  general  practice  to  calculate 
theoretical  thermodynamic  properties  for  a  proposed  formulation  based  on  its 
chemical  composition  and  the  heats  of  formation  of  its  ingredients.  In  this 
way , theoretical  values  for  impetus,  co-volume,  molecular  weight  of  combustion 
products  and  the  like  are  available  for  new  formulations  from  computations 
using  codes  such  as  BLAKE.  However,  these  theoretical  thermodynamic  values 
still  need  to  be  verified.  Fortunately,  while  impetus  is  strongly  sensitive 
to  composition,  co-volume  is  not.  As  stated  by  Hunt, 

"...,  since  impetus  is  sensitive  to  slight  variations 
in  propellant  composition,  particularly  in  the 
moisture  content,  the  experimental  values  of  impetus 
should  be  more  reliable  than  the  calculated  values. 

The  co-volume,  on  the  other  hand,  is  not  sensitive  to 
small  changes  in  composition  and  can  only  be 
derived  experimentally  by  a  difference  process;  the 
calculated  values  are  therefore  more  reliable." 

The  co -volumes  computed  by  BLAKE,  therefore,  can  be  assumed  as  reasonable  for 
most  systems,  except  for  those  having  high  levels  of  condensable  combustion 
products  (metal  oxides/  nitrides/  sulfides/  sulfates  and  their  derivatives). 
The  experimental  impetus  computations  described  below  are  predicated  on  the 
assumption  that  BLAKE  co-volumes  are  reasonably  correct  for  the  formulation. 

In  determining  experimental  impetus  a  sample  of  propellant  together  with 
an  appropriate  amount  of  igniter  is  completely  burned  in  a  closed  chamber  of 
known  volume.  The  pressure- time  history,  including  maximum  observed  pressure, 
is  recorded.  Based  upon  the  maximum  pressure,  the  closed  chamber  volume, 
thermochemical  properties  of  the  igniter,  and  co -volume  of  the  propellant  an 
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experimental  impetus  can  be  determined.  (Data  from  the  same  experiment  can 
also  be  used  for  burning  rate  computations).  The  method  works  best  for  fast¬ 
burning  samples  where  heat  loss  is  small. 

In  deriving  the  equation  for  computing  the  experimental  impetus  the 
following  assumptions  have  been  made  in  order  to  meet  the  objectives  stated  in 
the  introduction. 

1.  The  original  volume  of  the  closed  chamber  is  known. 

2.  Both  the  igniter  and  propellant  are  completely  burnt. 

3.  The  maximum  pressure  from  the  firing  is  known. 

4.  All  thermochemical  properties  of  the  igniter  are  known. 

5.  The  co-volume  of  the  propellant  is  known. 

6.  The  Noble-Able  equation  of  state  can  be  used  to  describe 
the  final  state  of  the  closed  chamber. 

7.  The  theoretical  impetus  of  the  igniter  can  be  used  in 
place  of  the  experimental  impetus  of  the  igniter. 

For  a  closed  chamber  with  free  volume,  Vb,  the  Noble-Able  equation  of 
state  is 


P(Vb  -  b)  =>  NRT 


(27) 


where  P  —  pressure 

b  =  co -volume  correction  term 
N  =  number  of  moles  of  gas 
T  -  system  temperature 
R  =  universal  gas  constant 

Since  all  igniter  and  propellant  is  assumed  to  be  completely  burnt,  the 
free  volume  of  the  chamber,  Vb,  is  the  volume  of  the  empty  chamber,  V.  Also, 
having  all  igniter  and  propellant  burnt  implies  that  the  co -volume  correction 
term,  b,  is  the  sum  of  the  product  of  the  co-volume  and  weight  for  the  igniter 
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and  propellant.  That  is, 


b  -  +  iipC  (28) 

Substituting  into  (27)  yields 

P(V  -  n^  -  n^C)  -  NRT  (29) 

Now  the  gas  in  the  chamber  is  a  mixture  of  gas  products  from  the  igniter  and 
the  propellant.  Thus,  N  -  +  Np  T:’.j.ch  when  substituted  into  (29)  produces 

the  following  equation. 

P(V  -  njCj  -  npC)  -  RT(N£  +  Np)  (30) 

Multiplying  the  right-hand  side  of  (30),  equation  (31)  is  obtained. 

P(V  -  n^  -  npC)  -  RTNi  +  RTNp  (31) 

Since  N  -  weight  divided  by  molecular  weight,  and  Np  can  be  replaced  to 
obtain 


P(V  -  n^  -  iipC)  -  RTC^/m^  +  RTC/nip  (32) 

Now  impetus  is  given  by  the  formula 

F  =  RT/'m  (33) 

Thus,  using  the  assumption  that  the  theorecical  impetus  of  the  igniter  can  be 
used  in  place  of  the  experimental  impetus  together  with  (33)  the  expression 
below  is  obtained. 


Finally,  equation  (34)  can  be  solved  for  the  experimental  impetus  of  the 
propellant. 


Fp  -  [P(V  -  n ^  -  tipC)  -  qF^/C 


(35) 


III.  MINICB  -  PROGRAM 

Program  Flow 

The  flow  of  the  program  MINICB  is  presented  below.  A  complete  listing 
of  the  program  together  with  a  listing  of  variables  used  in  it  can  be  found  in 
Appendix  B. 


END 
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Except  for  the  five  point  differentiation  the  calculations  used  in  MINICB 
were  described  in  section  II.  As  for  the  five  point  differentiation  method  it 
was  obtained  from  a  paper  by  Juhasz  and  Price ^  and  is  given  by  the  following 
formula. 

-2*P_2  -  P.i  +  0  +  P+1  +  2*P+2 

dPQ/dt  »  -  (36) 

10*  At 

In  the  formula  P_2>  P.^,  Pq,  P+^,  and  P+2  are  five  successive  points  in  the 
data  set  and  the  derivative  of  the  pressure  is  to  be  estimated  at  the  pressure 
value  of  Pq. 

RUNNING  MINICB 

To  run  MINICB,  load  the  program  MINICB  and  run  it  from  BASIC.  For  the 
most  part  MINICB  is  menu  driven  and  will  only  require  entering  the  necessary 
values  for  whatever  computation  is  to  be  performed.  The  input  values  needed 
for  each  option  are  listed  below  in  Tables  1  and  2  together  with  the  units  for 
each  input  value. 


TABLE  1.  Required  Input  for  the  Experimental  Impetus 
Calculation 

Flame  Temperature  of  Igniter  --  Degrees  Kelvin 

Molecular  Weight  of  Igniter 

Co -volume  of  Igniter  --  Cubic  Inch/lb. 

Weight  of  Igniter  --  Pounds 

Theoretical  Impetus  of  Igniter  --  FT. -lb. /lb.  Optional 

Weight  of  Propellant  --  Pounds 

Co-volume  of  Propellant  --  Cubic  Inch/lb. 

Observed  Maximum  Pressure  --  Psi 

Volume  of  Closed  Chamber  --  Cubic  Centimeters 


TABLE  2.  Required  Input  for  the  Burning  Rate 
Analysis 

Theoretical  Impetus  of  Igniter  --  Ft. -lb. /lb.  Optional 
Flame  Temperature  of  Igniter  --  Degrees  Kelvin 
Molcular  Weight  of  Igniter 
Co-volume  of  Igniter  --  Cubic  Inch/lb. 

Weight  of  Igniter  --  Pounds 

Weight  of  Propellant  --  Pounds 

Co-volume  of  Propellant  --  Cubic  Inch/lb. 

Flame  Temperature  of  Propellant  --  Degrees  Kelvin 

Molecular  Weight  of  Propellant 

Density  of  Propellant  --  lb. /Cubic  Inch 

Length  of  Propellant  Grain  --  Inch 
Diameter  of  Propellant  Grain  --  Inch 
Perf  Diameter  of  Propellant  Grain  --  Inch 
Inner  Web  of  Propellant  Grain  --  Inch 
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TABLE  2.  Required  Input  for  the  Burning  Rate 
Analysis  (Con't) 

Middle  Web  of  Propellant  Grain  --  Inch 
Outer  Web  of  Propellant  Grain  --  Inch 
Number  of  Perfs 
Form  Function 

Observed  Maximum  Pressure  --  Psi 

Volume  of  Closed  Chamber  --  Cubic  Centimeters 

Delta  Time  --  Msec 
Starting  Time  --  Msec 

Pressure  Data  --  Psi 


The  pressure  data,  at  fixed  time  intervals,  can  be  entered  by  hand  or 
read  from  a  file  stored  on  a  disk.  In  either  case,  the  pressure  values  are 
entered  one  per  line  with  the  last  value  being  a  -1  to  terminate  the  input. 
The  input  should  not  contain  time,  this  information  is  generated  from  knowing 
delta  time  and  the  starting  time.  Also,  the  pressure  should  always  be 
increasing.  If  the  pressure  is  not  increasing  then  the  program  will  produce 
negative  burning  rates  (physically  not  meaningful)  for  those  times  where  the 
pressure  decreases.  The  negative  bum  rates  reflect  the  fact  that  in  the 
analysis  a  decrease  in  pressure  means  an  increase  in  solid  propellant,  not 
heat  loss,  which  is  accounted  for  by  an  adjusted  system  temperature. 

Output  from  MINICB  is,  of  course,  dependent  upon  which  type  of 
calculation  is  being  performed.  For  the  experimental  impetus  calculation  the 
output  is  a  one  page  summary  of  the  input  values  together  with  the  computed 
experimental  impetus  for  the  propellant  given  in  Joules/gram  and  Ft. -lb. /lb. 
The  output  for  the  burn  rate  analysis  is  in  two  parts.  The  first  is  a 
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hardcopy  to  the  printer  which  consists  of  a  summary  of  the  input  values,  the 
theoretical  maximum  pressure,  number  of  propellant  grains  and  an  optional 
table  of  burning  rates  at  every  500  psi.  The  second  form  of  output  is  a  disk 
file  which  consists  of  a  table  containing  time,  pressure,  dP/dt,  percent 
burnt,  burning  rate,  and  surface  area  fraction.  This  file  will  also  contain 
the  table  of  burning  rates  at  every  500  psi  if  the  option  was  selected. 
Samples  of  the  output  available  from  MINICB  are  given  in  Appendix  C. 


IV.  VALIDATION 

Validation  of  MINICB  requires  an  assessment  of  its  performance  in 
computing  burning  rates  with  its  various  form  functions  as  well  as  in 
computing  experimental  impetus .  The  procedures  followed  and  results  obtained 
are  given  below. 


Burning  Rate  Validation 

O 

A  well- tested  lump  parameter  interior  ballistic  code,  IBHVG2,  was  used 
to  simulate  closed  bomb  burning  by  making  the  projectile  weight  so  large  that 
no  projectile  motion  occurred.  This  permitted  direct  quantitative  comparison 
of  MINICB  burning  rate  outputs  with  the  burning  rate  inputs  used  to  generate 
the  pressure-time  data  up  to  the  point  of  Pmax.  In  addition,  the  results 
from  MINICB  were  found  to  be  in  reasonably  close  agreement  with  experimental 
burning  rate  data  computed  by  CBRED^  (See  Appendix  G  and  Table  F-3,  Appendix 
F). 


At  the  present  time,  MINICB  contains  five  different  form  functions  for 
grain  geometry.  These  form  functions  are: 

1.  Cord  burning  cigarette  fashion  --  burning  on  one  end  only 

2.  Single  perf  burning  only  on  the  perf  surface 

3.  Single  perf  grain 
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4.  Generalized  n-th  perf  grain  --  the  form  function  assumes 
all  webs  are  equal  and  all  perf  diameters  are  equal.  In 
addition,  the  form  function  only  holds  up  to  the  point  of 
slivering. 

5.  Single  perf  burning  on  both  ends  and  the  perf  surface 

All  five  form  functions  were  exercised,  see  validation  runs  1-5  below.  For 
the  general  n-th  perf  grain  a  seven  perf  case  was  chosen.  More  details 
concerning  the  form  functions  can  be  found  in  Appendix  A. 

Propellants  A,  B,  and  C  with  bPn  burning  rate  laws  and  realistic 
thermochems  were  used  for  the  simulations.  The  values  of  b  and  n  as  well  as 
the  propellant/igniter  thermochemical  data  used  in  the  various  simulations  are 
given  in  Tables  5-5. 


TABLE  3.  Details  of  propellant  and  igniter  used 

in  IBHVG2  to  simulate  closed  bomb  firings 
for  validation  runs  1,  2  and  3. 


Propellant 


Igniter 


Type: 

A 

Coefficient  b 

0.0826 

Exponent  n 

0.865 

Weight (lb) : 

0.0165 

0.000001 

Density(lb/inA3) : 

0.056 

Gamma: 

1.2386 

1.2386 

Impetus (ft-lb/lb) : 

347243 

347243 

Co-volume(inA3/lb) : 

26.62 

26.62 

Flame  Temperature (K) : 

3092 

2000 
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TABLE  4.  Details  of  propellant  and  igniter  used  in  IBHVG2 
to  simulate  closed  bomb  firings  for  validation  4. 


Propellant 

Igniter 

Type: 

B 

Black  Powder 

Coefficient  b 

0.001216 

Exponent  n 

0.865 

Weight (lb) : 

0.1213 

0.00441 

Density(lb/inA3) : 

0.0596 

Gamma: 

1.2256 

1.25 

Impetus (ft- lb/lb) : 

361059 

96000 

Co -volume (inA 3/lb) : 

27.261 

30.00 

Flame  Temperature (K) : 

3289 

2000 

TABLE  5 .  Details 

of  propellant  and  igniter  used  in  IBHVG2 

to  simulate  closed  bomb  firings 

Propellant 

for  validation  5 

Igniter 

Type: 

C 

Black  Powder 

Coefficient  b 

0.001216 

Exponent  n 

0.865 

Weight(lb) : 

0.1213 

0.00441 

Density(lb/inA3) : 

0.0585 

Gamma: 

1.238 

1.25 

Impetus (ft- lb/lb) : 

359671 

96000 

Co -volume (inA3/lb) : 

28.482 

30.00 

Flame  Temperature (K) : 

3119 

2000 

Validation  1.  The  burning  of  propellant  A,  see  Table  3,  with  a  cord  geometry 
burning  cigarette  fashion  (i.e.  on  a  single  end  surface)  was  simulated  via 
IBHVG2  and  reduced  via  MINICB.  In  place  of  showing  all  the  output  from 
MINICB,  the  results  will  be  displayed  graphically  showing  the  percent 
difference  between  the  burning  rates  calculated  by  MINICB  and  the  actual 
burning  rates  as  given  by  the  burning  rate  law  r  -  .0826PA.865  for  different 
values  of  pressure,  see  Figure  2.  Agreement  between  the  burning  rate  values 
extracted  and  the  initial  burning  rate  values  was  close  for  the  entire  burn, 
with  a  mean  error  of  0.57  percent. 

CIGARETTE  FORM  FUNCTION 


PBBCB5T  DIPPEBE5C2 


Figure  2 .  Percent  difference  between  computed 

burn  rates  from  MINICB  and  theoretical 
burning  rates  for  a  cord  burning  in 
a  cigarette  fashion.  Validation  1. 
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PERCENT  DIFFERENCE 


Validation  2:  Propellant  A,  Table  3,  single  perforated  grain,  burning  only 
on  the  perforation  surface.  Results  similar  to  the  first  case  are  presented 
below  in  Figure  3.  In  this  case  the  mean  error  was  -.87  percent. 


BURNING  ON  PERF  SURFACE  ONLY 


PBBCBRT  DIFFBRBHCB 


Figure  3.  Percent  difference  in  burning  rates  between 
MINICB  and  rates  computed  using  the  burning 
rate  law  for  a  single  nerf  grain  burning  on 
the  nerf  surface  only,  validation  2. 
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Validation  3.  Propellant  A,  single  perforated  geometry,  burning  on  the  ends 
and  the  perforation  surface.  Results  similar  to  those  above  are  given  in 
Figure  4.  In  this  case  the  mean  error  was  1.69  percent.  The  error,  however, 
was  increasing  near  the  end  of  the  burn,  reaching  a  maximum  of  4  percent  at 
burnout . 


BURNING  ON  PERF  6c  END  SURFACE 


PBBCKKT  DIFFBBBIfCB 


Figure  4 .  Percent  difference  in  burning  rates  between 
MINICB  and  rates  computed  using  the  burning 
rate  law  for  a  single  perf  grain  burning  on 
the  ends  and  oerf  surface,  validation  3. 


Validation  5.  Propellant  C,  Table  5,  seven  perforated  grain  burning  on  all 
surfaces.  Results  are  presented  in  the  same  manner  as  above  in  Figure  6.  The 
burning  rates  computed  were  approximately  two  percent  high  initially  and  1.75 
percent  low  at  the  end  of  the  analysis  which  was  performed  to  the  point  of 
grain  slivering.  The  overall  absolute  mean  error  was  0.88  percent. 


SEVEN  PERF  GRAIN 


Figure  6 .  Percent  difference  in  burning  rates  between 
MINICB  and  rates  computed  using  the  burning 
rate  law  for  a  seven  oerf  grain,  validation  5. 


Overall,  for  the  different  form  functions,  the  average  error  is  less  than 
two  tentns  of  one  percent  in  the  best  case  and  below  1.4  percent  in  the  worst 
case.  On  the  5  to  95%  interval  of  Pmax  the  greatest  observed  error  was  -3%. 

Due  to  the  lower  level  of  error  in  the  constant  area  case ,  it  seems  probable 
that  numerical  errors  are  magnified  during  the  surface  area  computations .  The 
error  levels  above  fall  within  the  normal  experimental  round- to -round 
variation  of  closed  bomb  burning  rates  and,  therefore,  are  not  overly 
critical.  Considering  the  simplifying  assumptions  made  in  developing  the 
theory  and  the  crude  numerical  methods  used,  the  error  levels  appear 
reasonable  for  the  intended  application.  MINICB  was  designed  as  a  quick  tool 
for  determining  burning  rates  from  a  limited  number  of  pressure  points.  The 
simplifications  introduced  were  necessary  to  permit  running  in  a  minimum 
environment. 


Burn  Rate  Error  vs.  Number  of  Data  Points 

Since  MINICB  can  be  run  using  a  programmable  calculator  or  a  computer 
with  pressure  information  entered  by  hand,  it  was  of  interest  to  investigate 
what  effect  the  number  of  pressure  values  entered  would  have  on  program 
performance.  The  data  set  for  the  end-burning  cord,  see  verification  run  1 
for  details,  was  reduced  via  MINICB  using  different  numbers  of  pressure 
values.  In  all  instances,  points  were  selected  at  equal  intervals  and  included 
the  original  first  pressure  value  as  well  as  the  maximum  pressure  value.  Data 
sets  containing  253,  126,  64,  22,  or  12  points  were  used.  The  results  are 
presented  in  terms  of  percent  difference  and  mean  and  standard  deviation  of 
percent  difference,  see  Figure  7  and  Table  6. 
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PEBCENT  DIPPEBEMCB 


CIGARETTE  FORM  FUNCTION 


PERCBNT  DIFFERENCE 


A  22  POINTS  X  18  POINTS 

Figure  7 .  Percent  Difference  in  burning  rates  between 

MINICB  and  theoretical  burning  rates  for  a  cord 
burning  in  a  cigarette  fashion  for  data  sets 
with  different  number  of  points  where  each 
set  was  taken  from  the  same  original  data  set. 


TABLE  6.  Mean  and  standard  deviation  of  percent 

error  for  data  sets  with  different  number 
of  points  all  of  which  were  taken  from 
the  same  original  data  set. 

Number  of  Points  Mean  Standard  Deviation 


253 

126 

64 

22 

12 


.567233 

.474999 

.504159 

.325364 

.42504 

.093079 

.274162 

.074713 

.17125 

.086029 
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Paradoxically,  as  the  number  of  data  points  in  the  set  decreased  the  mean 
percent  error  decreased.  This  is  probably  related  to  numerical  roundoff  errors 
in  the  calculated  mass  of  propellant  burnt  for  the  smaller  time  steps. 

Similar  examinations  using  the  other  form  functions  showed  little  or  no 
dependence  of  the  percent  error  on  the  number  of  data  points  utilized. 
Ultimately,  this  indicates  that  the  method  should  be  reasonably  reliable  when 
using  a  small  number  of  manually  entered  pressure  points  to  compute  burn 
rates . 


Experimental  Impetus 

It  was  not  practical  to  validate  experimental  impetus  calculations  via 
test  firings.  It  was  possible,  however,  to  verify  the  calculation  by 
comparison  with  results  from  the  BLAKE  thermodynamic  code.  BLAKE  computes 
pressure,  impetus,  co-volume,  molecular  weight,  and  other  thermodynamic 
properties  of  propellant  combustion  products  at  various  loading  densities.  In 
essence,  the  pressure  from  a  BLAKE  computation  at  a  given  loading  density  is 
the  equivalent  of  the  maximum  pressure  from  an  adiabatic  closed  bomb  test. 
Under  the  circumstances,  the  experimental  impetus  computed  by  the  proposed 
method  should  be  the  same  as  the  theoretical  impetus  from  BLAKE. 

The  propellant  composition  chosen  for  the  evaluation  was  a  ball 
propellant,  X-4179,  manufactured  by  the  Olin  Corporation.  The  composition  of 
X-4179  is  given  in  Figure  8.  The  input  to  BLAKE  is  given  in  Figure  9  and  the 
output  from  BLAKE  appears  in  Figure  10. 
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To  Whom  It  May  Concern: 


Shown  below  are  the  analytical  results  on  the  following  unrolled  pro¬ 
pellants: 


X-4179 

X-4180 

X-4182 

X-4183 

Coating  Series  No. 

340-104 

223-33 

418-143 

208-193 

Ave.  Grain  Diameter 

.0284“ 

.0238“ 

.0271“ 

.0163" 

Nitrogen  of  NC,  X 

13.16 

13.02 

13.04 

13.16 

Nitrocellulose,  X 

59.94 

65.27 

75.17 

85.95 

Nitroglycerin,  X 

37.86 

32.54 

22.77 

11.91 

Dibutylphthalate,  X 

0.28 

0.25 

0.24 

0.37 

Diphenylamine,  X 

0.49 

0.84 

0.84 

0.77 

N-Nitroso  Diphenylamine,  X 

0.81 

0.53 

0.53 

0.54 

Dinitrotolvene,  X 

0.14 

0.14 

0.08 

0.07 

Water,  X 

0.46 

0.39 

0.44 

0.77 

Moisture  and  Volatiles,  X 

0.52 

0.43 

0.51 

0.94 

Total  Volatiles,  X 

0.51 

0.42 

0.50 

1.01 

Residual  Solvent,  X 

0.04 

0.03 

0.06 

0.24 

Ethyl  Acetate,  X 

0.00 

0.00 

0.03 

0.22 

Benzene,  X 

0.04 

0.03 

0.03 

0.02 

Sodium  Sulfate,  X 

0.09 

0.09 

0.08 

0.09 

Calcium  Carbonate,  X 

0.04 

0.04 

0.03 

0.02 

Ash,  X 

0.15 

0.13 

0.09 

0.17 

Graphite,  X 

0.20 

0.17 

0.17 

0.11 

Dust,  X 

0.03 

0.04 

0.00 

0.00 

120°  Stability,  SP  Time,  Min. 

55 

70 

80 

115 

*idan 

(Technical  Services 
peration 


DZ6.014.2 


Figure  8 .  Composition  of  X-4179.  Olin  Corporation 
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X-4179 


THE  COHPOSETION  IS 


NAHE 

PCT  NT 

PCT  MOLE 

DEL  H-CAL/H 

FORMULA 

NC1316 

58,986 

.087 

-1.6457EI0B 

C 

H 

0 

N 

6000 

7361  10278 

2639 

NB 

37.258 

67.996 

-8.8600E*04 

C 

H 

0 

N 

3 

5 

9 

3 

OBP 

.276 

.411 

-2. 0140E+05 

C 

H 

0 

16 

22 

4 

OPA 

.433 

1.182 

3.1070E<04 

C 

N 

H 

12 

1 

11 

SLOP 

.789 

1.649 

5.0930E«04 

C 

H 

N 

0 

12 

10 

2 

1 

ONI 

.138 

.314 

-1.7100E*04 

C 

K 

0 

N 

7 

6 

4 

2 

H20 

.460 

10.576 

-6.8315E*04 

H 

0 

2 

1 

ETOH 

1.069 

9.620 

-6.6420Et04 

C 

H 

0 

2 

6 

1 

C6H6 

.039 

.207 

1.1720E*04 

C 

H 

6 

6 

SODStl 

.089 

.260 

-2.6920E«05 

NA 

S 

0 

2 

I 

4 

CALC 

.217 

.907 

-2.8851E*05 

K 

C 

0 

I  1  3 


C  .197  6.793  0.  C 

I 

THE  HEAT  DF  FORMATION  IS  -530.95  CAL/GRAH  »  -2.2005E+05  CAL /HOLE. 

THE  ELEMENTS  AND  PERCENT  BY  MOLE 
C  20.577 

H  23.426 

0  39.626 

N  11.327 

NA  .013 

S  .007 

K  .023 


Figure  9.  Composition  for  X-4179  used  in  BLAKE 
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mil  TEMP  PRESSURE  IMPETUS  HOL  NT  CO-VOL  FROZEN  CP(FR)  BIT)  C1T)  GftS  VOL  S  H  E  ADEXP  PHI 
6/CC  K  PSI  FT-LB/LB  6fl5  CU  IN  BWWfl  CAUH-K  CU  IN  IN«»6  CD  IN/LB  6IDBS  CAL/B  CfiL/G 


1) 

.0500 

3534. 

8743. 

362444. 

24.049 

28.38 

1.2145 

11.30 

2) 

.1000 

3435. 

18424. 

384443. 

24.142 

27.83 

1.2139 

11.39 

J) 

.1500 

3441. 

29447. 

388500. 

24.214 

27.28 

1.2144 

11.45 

4) 

.2000 

3460. 

41898. 

389901. 

24.250 

26.73 

1.2155 

11.52 

5) 

.2500 

3493. 

55479. 

390944. 

24.277 

26.15 

1.2172 

11.58 

6) 

.2000 

3704. 

70494. 

391813. 

24.299 

25.57 

1.2196 

11.64 

7) 

.2500 

3714. 

87047. 

392515. 

24.317 

24.97 

1.2225 

11.71 

8) 

.4000 

3722. 

105257. 

393113. 

24.333 

24.38 

1.2260 

11.78 

9) 

.4500 

3729. 

123204. 

393427. 

24.348 

23.78 

1.2299 

11.85 

10) 

.5000 

3735. 

147014. 

394073. 

26.341 

23.19 

1.2344 

11.93 

11) 

.5500 

3740. 

170787. 

394442. 

26.374 

22.61 

1.2393 

12.02 

12) 

.4000 

3745. 

194431. 

394803. 

26.386 

22.04 

1.2446 

12.10 

13) 

.4500 

3750. 

224651. 

395100. 

26.399 

21.48 

1.2504 

12.19 

HI 

.7000 

3754. 

254952. 

395357. 

26.411 

20.93 

1.2565 

12.2? 

15) 

.7500 

3758. 

287439. 

395578. 

24.424 

20.40 

1.2630 

12.39 

14) 

.8000 

3741. 

322814. 

295762. 

26.437 

19.89 

1.2699 

12.49 

1.652 

2.19 

553.598 

2.29 

-242.8 

-531.0 

1.2399 

I.054U 

1.650 

2.19 

276.799 

2.23 

-224.0 

-530.9 

1.3141 

1.1118 

1.648 

2.19 

184.533 

2.20 

-205.2 

-531.0 

1.3871 

1.1735 

1.647 

2.19 

138.399 

2.17 

-185.7 

-531.0 

1.4601 

1.2393 

1.647 

2.19 

110.720 

2.15 

-165.3 

-531.0 

1.5331 

1.3093 

1.646 

2.19 

92.266 

2.14 

-143.7 

-531.0 

1.6059 

1.3833 

1.645 

2.19 

79.085 

2.12 

-121.1 

-531.0 

1.6782 

1.4615 

1.645 

2.19 

69.195 

2.11 

-97.4 

-531.0 

1.7498 

1.5439 

1.644 

2.19 

61.508 

2.10 

-72.5 

-531.0 

1.8205 

1.6304 

1.644 

2.19 

55.358 

2.09 

-46.4 

-531.0 

1.8903 

1.7210 

1.644 

2.19 

50.326 

2.08 

-19.3 

-531.0 

1.9590 

1.8158 

1.643 

2.19 

46.132 

2.08 

9.1 

-531.0 

2.0266 

1.9147 

1.643 

2.18 

42.584 

2.07 

38.6 

-531.0 

2.0931 

2.0177 

1.643 

2.18 

39.542 

2.06 

69.2 

-531.0 

2.1585 

2.1250 

1.643 

2.18 

36.906 

2.06 

101.0 

-531.0 

2.2228 

2.2363 

1.643 

2.18 

34.600 

2.05 

134.0 

-531.0 

2.2859 

2.351B 

Figure  10 .  Results  from  BLAKE  for  X-4179 


Using  the  BLAKE  results  for  a  loading  density  of  0.25  g/cc  with  a  bomb 
volume  of  100  cc  the  weight  of  the  igniter  (assumed  to  have  identical 
thermochemical  properties  to  the  propellant)  is  selected  as  1  gm  and  that  of 
the  propellant  24  gm.  Figure  11  shows  the  results  of  MINICB  using  the 
experimental  impetus  option  for  this  propellant  with  input  values  for  the 
necessary  thermochemical  values  taken  from  BLAKE. 
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IMPETUS  CALCULATION 


IGNITER  IMPETUS: 
IGNITER  CO -VOLUME: 
IGNITER  WEIGHT: 


390964  FT-LB/LB 
26.15  CU  IN/LB 
.002205  LB 


PROPELLANT  WEIGHT:  .05292  LB 

PROPELLANT  CO- VOLUME:  26.15  CU  IN/LB 


BOMB  VOLUME: 
MAXIMUM  PRESSURE: 


100  CC 
55479  PSI 


IMPETUS-  1163. 414  J/G 
IMPETUS-  390895  FT-LB/LB 


Figure  11.  Experimental  Impetus  results  from 

Ml_NICB..£or  .X-.4.179. 


From  Figure  10  the  theoretical  impetus  for  a  loading  density  of  0.25  g/cc 
is  390964  Ft-lb/lb.  MINICB  calculates  a  impetus  of  390895  Ft-lb/lb  which 
represents  a  percent  error  of  -0.0176%.  This  small  difference  is  probably  a 
result  of  round-off  errors.  Considering  the  purpose  of  the  computations,  which 
is  to  compute  experimental  impetus  as  a  practical  performance  index,  the 
agreement  is  reasonable. 


V.  SUMMAPvY 

MINICB,  with  overall  error  limits  of  at  most  1.5%,  can  be  a  particularly 
useful  tool  ir  the  computation  of  burning  rates  or  experimental  Impetus  from 
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closed  bomb  pressure- time  curves.  This  is  especially  true  in  situations  where 
only  a  limited  number  of  pressure- time  data  points  are  available,  the 
propellant  has  an  unusual  grain  geometry,  or  is  a  new  formulation  and  feedback 
is  quickly  needed  for  the  experimental  program  to  proceed.  Considering  the 
minimal  environment  necessary  for  running  MINICB  and  the  simplicity  of  use,  it 
is  felt  that  MINICB  meets  the  stated  objects  set  out  in  the  introduction. 
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TA3LE  OF  SYMBOLS 


L 


“p 


"p 


n„- 


Pmax 

Ptheo 


Instantaneous  propellant 
surface  area 

Co -volume  correction  term 

Weight  of  the  Propellant 

Weight  of  Propellant  Burnt 

Weight  of  the  Igniter 

Weight  of  Igniter  Burnt 

Perf  diameter 

Propellant  outer  diameter 

Experimental  Impetus  of 
propellant 

Theoretical  impetus  of  Igniter 
Grain  length 

Average  molecular  weight  of 
Propellant  combustion  products 

Average  molecular  weight  of 
Igniter  combustion  products 

Number  of  moles  of  gas 

Number  of  moles  of  igniter 
gas  products 

Number  of  moles  of  propellant 
gas  products 

Co -volume  of  the  Propellant 
Combustion  products 

Co-volume  of  the  Igniter 
Combustion  products 

Observed  maximum  pressure 

Theoretical  maximum  pressure 
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TABLE  OF  SYMBOLS  (Con't) 


Pi 


Pp.P 

Pi 

P 


R 


r 

T 


T 


i 


x 

V 


Mathematical  pi 
Density  of  the  Propellant 
Solid  igniter  density 
Chamber  Pressure 
Universal  Gas  Constant 
Burn  rate 

System  gas  temperature 

Flame  Temperature  of  the 
Propellant 

Flame  Temperature  of  the 
Igniter 

Linear  distance  burnt 
Volume  of  empty  chamber 
Free  Volume  in  the  bomb 
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APPENDIX  A 

SURFACE  AREA  FORMULAS 
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I.  Cord  Burning  Cigarette  Fashion:  An  example  of  this  type  of  burning  grain 
is  shown  in  Figure  A-l. 


Figure  A-l.  A  cord  burning  in  a  cigarette  fashion. 


For  this  particular  grain  geometry  the  burning  is  all  on  one  end  and  the 
surface  area  of  the  burning  surface  remains  constant  and  is  given  by  equation 
(A-l). 


A  =  (Pi/4)*D2 


(A-l) 


II.  Single  Perf  Burning  Only  on  Perf  Surface:  A  diagram  of  the  way  in  which 
the  grain  is  burning  is  shown  in  Figure  A- 2. 
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The  equation  for  the  volume  burnt  is  given  in  equation  (A-2) 


Volume  Used  -  L*Pi*[((d  +  2x)/2)2  -  (d/2)2] 


Squaring  and  simplifying  in  the  brackets  yields 


Volume  Used  =  L*Pi*[dx  +  x2] 


Thus,  multiplying  by  the  density,  which  is  denoted  by  p  instead  of  Pp,  to 
obtain  the  total  propellant  burnt,  equation(A-4)  is  obtained. 


Cp  -  p*L*Pi*[dx  +  x2] 


(A-4) 


Rearranging,  the  following  quadratic  equation  in  x  is  obtained. 


x2  +  dx  -  Cp/(p*Pi*L)  -  0 


(A-5) 


Letting  B  -  Gp/(p*Pi*L)  and  using  the  quadratic  formula  to  solve  for  x. 


x 


d 


+  4*B 


2 


(A-6) 


Since  x  is  positive  the  positive  root  is  chosen  in  the  expression  for  x.  Next 
the  equation  for  the  surface  area  for  this  grain  is 


A  L*Pi*[d  +  2x] 


(A- 7) 


Substituting  (A-6)  into  (A-7)  yields  the  final  expression  for  the  surface  area 


A  -  L*?i*-y[<r 


+  4*Cp/(L*Pi*p) ] 
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(A-8) 


III.  Single  Perf  Burning  on  Ends  and  Perf  Surface:  A  diagram  of  this  grain  is 
shown  ir.  Figure  A- 3. 


Figure  A- 3.  Single  oerf  burning  on  ends 
and  perf  surface. 


Again  start  with  an  expression  for  the  total  volume  used. 


Volume  Used  =  L*Pi*[((d  *  2x)/2)2  -  (d/2)?i 

+  2x*Pl*[  (D/2)2  -  ((d  -i-  2x)/2)2]  (A-9) 


In  equation  (A-9)  the  first  term  represents  the  volume  burnt  from  inside  the 
perf  and  the  second  term  that  burnt  off  both  ends  which  was  rot  accounted  for 
in  the  first  term.  Simplifying  and  multiplying  by  density  to  obtain  the  total 
propellant  burnt  equation  (A-10)  is  obtained. 
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cp  -  (Pi*p/2)*[2*L*d*x  +  2*L*x2  +  xD2 
-  xd2  -  4*d*x2  -  4*x3] 


(A-10) 


Rearranging  to  obtain  an  equation  in  terms  of  x, 


x3  +  (-L/2  +  d)*x2  +  ( - L*d/2  -  D2/4  +  d2/4)*x 

+  Cp/(2*Pi*p)  -  0  (A- 11) 


Now  this  is  a  cubic  equation  in  x  and  the  tendency  is  to  assume  that  there  is 
one  real  root  and  two  complex  roots.  However,  this  is  not  the  case.  To  see 
this,  note  that  the  intercept  is  positive  and  that  the  lead  coefficient  also 
positive.  Thus,  if  a  graph  was  made  of  this  function  it  would  tend  to 
positive  infinity  as  x  went  to  positive  infinity  and  would  be  positive  for  x 
equal  to  zero,  see  Figure  A-4.  Therefore,  the  graph -must  cross  the  positive 
x-axis  twice  or  not  at  all.  However,  if  it  does  not  cross  the  positive  x-axis 
then  no  solution  would  exist,  but  one  must  exist.  Hence,  equation  (A-ll)  has 
three  real  roots,  two  positive  and  one  negative. 
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Figure  A- 4.  Graph  of  the  Equation  Defined 
by  Equation  (A-ll) 


After  analyzing  solutions  of  cubic  equations  it  became  obvious  that  a 
numerical  solution  for  equation  (A- 11)  was  the  only  feasible  method  for 
solving  for  x  and  keeping  to  the  objective  of  having  a  simple  and  easy  to  run 
program.  The  Newton  Iterative  method,  starting  with  an  initial  guess  of  0, 
for  finding  roots  has  proven  sufficient  in  computing  the  root. 

Having  obtained  a  value  for  x,  the  surface  area  is  giver,  by  equation  (A-12),, 


A  =  Pi*[d  +  2x]  +  ?.*Pi*[(D/2)2  -  ((d  +  2x)/2)2] 


(A-12) 
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IV.  N-th  Perf  Grain:  An  example  of  this  grain  geometry  for  seven  perfs  is 
shown  in  Figure  A- 5. 


(Typ.  7  places) 
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Figure  A- 5.  S^ven  Pevf  Grain 


In  the  following  derivation  for  the  surface  area  the  analysis  is  only 
correct  up  to  the  point  of  slivering  and  assumes  that  the  webs  are  equal  and 
that  all  perf  diameters  are  the  same.  A  equations  will  be  derived  for  n  perfs 
not  just  rhe  seven  shown  in  Figure  A-5.  Also,  this  analysis  will  not  be  valid 
for  a  single  perf  grain.  As  before,  the  starting  point  is  the  volume  used  and 
is  given  in  equation  (A- 13)  which  is  very  similar  to  that  for  a  single  perf 
grain. 
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Volume  Used  -  L*Pi*[(D/2)2  '  2x)/2>  3 

+  n*L*Pi*[ ( (d  +  2x)/2)2  -  (d/2)2] 

+  2x*Pi*[ ( (D  -  2x)/2)2  -  n*( (d  +  2x)/2)2]  (A- 13) 


Multiplying  by  density  to  obtain  the  total  amount  of  propellant  burnt  and 
rearranging  to  obtain  an  equation  in  x,  equation  ^A-14)  is  obtained. 


(4n  -  4)x3  +  (4D  +  4nd  -  2Ln  +  2L)x2 

+  (nd2  -  D2  -  2Ldn  *  21.3 )x  +  2Cp/(Pi*p)  -  0 


(A-14) 


As  in  III.  the  cubic  also  has  three  real  roots  and  a  numerical  method  should 
be  used  to  determine  the  value  of  x.  Having  determined  x  the  surface  area  is 
given  by  equation  (A- 15). 


A  -  (Pi/2)*[(D  -  2x) 2  -  n(d  +  2x>2] 


+  Pi*(L  -  2x)*[n*(d  +  2x)  (D  -  2x)  ]  <A-15) 
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APPENDIX  B 

PROGRAM  LISTING  FOR  MINICB 
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PROGRAM  LISTING 
IBM  BASIC  VERSION 


10  REM  data  of  version:  17  december  1986 


100 

110 

120 

130 

140 

160 

160 

170 

180 

130 

200 

21° 

220 

220 


WIDTH  "LPT1 : " , CO : OPEN  "LPTi : "  FOR  OUTPUT  AS  #3 
DIM  R( 1028 , 9) 

COLOR  14 

REM  mini  closed  bomb  program  written  by  william  oberle 
REM  **  main  menu  and  selection  of  options  ** 

REM  ***************************v********************** 
COS: PRINT  : PRINT  : PRINT  : PRINT  : PRINT  : PRINT  PRINT 
PRINT  "  Mini  Closed  Bomb  Program" 

PRINT  "  OPTIONS" 

PRINT 

PRINT  "  1.  Perform  impetus  calculations" 

PRINT  "  2.  Reduce  pressure/time  data  to  burning  races" 

PRINT  M  3 .  End  this  program" 

PRINT  : PRINT 


2^0  INPUT  'Select  option  by  number" ;0 
2f-0  IF  0-3  THEN  END 

260  IF  (Ool  AND  002)  THEN  BEEP:  PRINT  "Not  a  selection" : GOTO  240 
270  GOSUB  190 : IF  0=1  THEN  GOSUB  740  ELSE  GOSUB  930 
280  GOTO  160 

290  REM  **  subroutine  to  enter  needed  information  for  options  1  and  2  ** 

300  REM  ****************************************************************** 
310  FT**0 :  CLS :  PRINT  :  PRINT  :  PRINT  :  PRINT  :  FRINT  :  PRINT 

320  PRINT  "Certain  information  is  needed  to  perform  the  selected  option." 
330  PRINT  "You  will  be  asked  for  this  information,  please  be  sure  to  enter" 
340  PRINT  "the  information  with  the  correct  units.  DO  NOT  HIT  RETURN!!" 

350  FOR  KK-1  TO  8000: NEXT  KK 
360  CLS 

370  PRINT  : PRINT  : PRINT  .’PRINT  : PRINT  : PRINT  : PRINT 
380  PRINT  "  IGNITER  INFORMATION" 

390  PRINT  : PRINT  "IF  THE  VALUE  OF  THE  IMPETUS  OF  THE  IGNITER  IS  KNOWN, ENTER 
ITS'' 


400  FRINT  "VALUE  IN  FT. -LB. /LB.  IF  BLACK  POWDER  IS  BEING  USED  ENTER  -1," 
410  PRINT  "OTHERWISE  HIT  RETURN." 

420  INPUT  FI:FI-FI*12 

430  IF  FI— 12  THEN  FI=97500! :TI-2170:FI-FI*12:M1-33.64;NI-21.8:G0T0  500 
440  IF  (0-1  AND  FIOO)  THEN  4  90 

*50  INPUT  "ENTER  THE  FLAME  TEKPSRATURF  OF  THE  IGNITER (K) " ,TI 
460  INPUT  "ENTER  THE  MOLECULAR  WEIGHT  OF  THE  IGNITE?." ; MI 
470  IF  FIOO  THEN  490 

430  f 1-8 . 31433*TI/MI/ . 0029890706# :  IT=FI*12 

490  INPUT  "ENTER  THE  CO- VOLUME (CU  IN/LB)  OF  THE  IGNITER" ;NI 

500  INPUT  "ENTER  THE  WEIGHT (LB)  OF  THE  ICN1TER";CI 

510  CLS 

520  PRINT  : PRINT  : PRINT  ’PRINT  : PRINT  : PRINT 
530  PRINT  "  PROPELLANT  DATA" 
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540  INPUT  "ENTER  THE  WEIGHT(LB)  OF  THE  PROPELLANT" ; CP 

550  INPUT  "ENTER  THE  CO-VOLUME(CU  IN/LB)  OF  THE  PROPELLANT" ;NP 

550  IF  0-1  THEN  680 

570  INPUT  "ENTER  THE  FIAME  TFMPERATURE(K)  OF  THE  PROPELLANT" ; TP 
580  INPUT  "ENTER  THE  MOLECULAR  WEIGHT  OF  THE  PROPELLANT" ;MP 
590  INPUT  "ENTER  THE  DENSITY  OF  THE  PROPELLANT (LB/CU  IN)";RHO 
600  CLS 

610  PRINT  • PRINT  : PRINT  : PRINT  : PRINT  : PRINT  : PRINT 
620  PRINT  "  GRAIN  INFORMATION" 

630  INPUT  "ENTER  THE  GRAIN  DIAMETER(IN) " ;D 

640  INPUT  "ENTER  THE  GRAIN  LENGTH (IN) " ;L 

650  INPUT  "ENTER  THE  NUMBER  OF  PERFS";NNP 

660  INPUT  "ENTER  THE  AVERAGE  PERF  DIAMETER (IN) " ; PD 

670  INPUT  "ENTER  THE  OUTER,  MIDDLE,  AND  INNER  WEB(IN) " ; WO, WM.WI 

680  CLS 

690  PRINT  : PRINT  : PRINT  -.PRINT  : PRINT  : PRINT  : PRINT 
700  PRINT  BOMB  AND  PRESSURE  INFORMATION" 

710  INPUT  "ENTER  THE  OBSERVED  PMAX(PSI) '* ;  PMAX 
720  INPUT  "ENTER  THE  BOMB  VOLUME (CC) " ;VB 
730  RETURN 

740  REM  **  IMPETUS  CALCULATIONS  ** 

750  REM  ****************************************************************** 
760  LPRINT  "  IMPETUS  CALCULATION" 

770  LFRINT  : LPRINT  : 

780  LPP.INT  "IGNITER  IMPETUS:  "; FI/12:"  FT-LB/LB" 

790  LPRINT  "IGNITER  CO- VOLUME:  ";NI;"  CU  IN/LB" 

800  LPRINT  "IGNITER  WEIGHT:  ";CI;"  LB" 

810  LFRINT 

820  LPRINT  "PROPELLANT  WEIGHT:  ";CP;M  LB" 

830  LPRINT  "PROFELIANT  CO-VOLUME:  ";NP;"  CU  IN/LB" 

840  LPRINT 

850  LPRINT  "BOMB  VOLUME:  ";VB;"  CC" 

860  LPRINT  "MAXIMUM  PRESSURE:  " ;PMAX; "  PSI" 

870  LPRINT  : LPRINT 
880  FI-FI/12 

890  FP-( (PMAX/145 .038)*(VB-CP*16 . 387045#*NP-CI*16 . 387045#*NI) 
-FI*.0029890706#*CI*453 . 592)/(CP*453 . 592) 

900  LPRINT  "IMPETUS-" ;FP; "  J/G" 

910  LPRINT  "IMPETUS-" ;INT(FP/. 0029890706#) FT-LB/LB" 

920  RETURN 

930  REM  **  BURNING  RATE  ANALYSIS  ** 

940  REM  *****************************************'*'****,*'************'**'!V***'**x 
950  PTHEO- ( 33372 ! *TP*CP/MP+FI*CI ) / (VB/16 . 38706  - CP*NP- CI*NI ) 

960  TOS-(CP*TP+CI*TI)/(CP+CI) 

970  TS— TOS*PMAX/PTHEO 

980  CONI— VB/16 . 38706 -CP/RHO-CI*NI 

990  C0N2-CI*FI*TS/TI 

1000  C0N3-33372 !*TS/MP 

1010  C0N4-NP - 1 /RHO 

1020  PIG=FI*TS*CI/TI/(VB/16.38706-CP/RH0-CI*NI) 
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1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1X2G 

1130 

U40 

1150 

1160 

11/0 

1130 

1190 

1200 

1210 

1220 

123C 

1240 

1250 

1260 

1270 

1200 

1290 

1300 


CLS: PRINT  : PRINT  : PRINT  : PRINT  : PRINT  : PRINT 

PRINT  "THIS  PROGRAM  ASSUMES  THAT  THF.  IGNITER  IS  ALL  BURNT  BEFORE  THE  " 
PRINT  "PROPELLANT  STARTS  TO  BURN.  THUS,  THE  POINTS  TO  BE  ENTERED  IN" 
PRINT  "THE  PRESSURE  VS  TIME  TABLE  MUST  START  UITH  A  PRESSURE  GREATER" 
PRINT  "THAN  THE  PRESSURE  DUE  TO  THE  IGNITER.  THIS  PRESSURE  IS  ";PIG;" 
PS  I" 

PRINT  : PRINT 

PRINT  "THE  PROGRAM  ALSO  ASSUMES  THAT  THE  TIME  STEP  IS  A  CONSTANT." 
INPUT  "PLEASE  ENTER  THIS  CONSTANT  TIME  STEP (MSEC) " ;DT:DT«DT/1000 
INPUT  "ENTER  THE  STARTING  TIME(MSEC) " ;TSTART:TSTART»TSTART/1000 
CLS 

PRINT  : PR?  NT : PRINT  : PRINT  : PRINT  : PRiNT 

PRINT  "THIS  PROGRAM  WILL  SUPPORT  THE  FOLLOWING  FORM  FUNCTIONS." 

PRINT  :  PRINT 


PRINT  " 

1. 

CIGARETTE" 

PRINT  " 

2, 

INHIBITED  TINGLE  PEF.F( BURNS  ONLY  ON  PLRF  SURFACE)" 

PRiNT  " 

3, 

SINGLE  PERF" 

PRINT  " 

4 

MULTI -PERF" 

PRINT  " 

5.. 

SINGLE  PLRF  INHIBITED  ON  LATERAL  OUTER  SURFACE" 

PRINT  : PRINT 

INPUT  "ENTER  YCUR  CHOICE  TOR  THE  FCRM  FUNCTION  BY  NUMBER" ;FFO 
CLS: PRINT  : PRINT  : PRINT  : PRINT  : PRINT  : PRINT 

PRINT  "THE  PRESSURE  DATA  MAT  BE  ENTERED  FROM  A  DATA  FILE  OR  BY  HAND." 
PRINT  "IF  A  DATA  FILL  IT  TO  BE  USED  THE  LAST  PRESSURE  MUST  BE  -X." 
PRINT  "THE  FILE  IS  l’C  TON"Y  fN  PRESSURE  VALUES  ONLY." 

PRINT  : PRINT  ‘INPUT  " ENTFR  T  FOR  FILE  INPUT  OR  H  FOR  HAND  INPUT" ;Z$ 
IT  ;i$="H"  J.GiN  1400 

IF  Z$0"F’-  THEN  PRINT  "NOT  A  CHOICE": GOTO  1240 
PRINT: PRINT 


1310  PRINT  "ENTER  IBP  NAME  OF  THL  FILE  CONTAINING  THE  PRESSURE  DATA." 
1320  PRINT  "2E  SURF.  TO  INCLUDE  THE  DRIVE  SPECIFICATION  IF  IT  IS  NOT  C." 
1330  INPUT  "THE  FILE  NAME  IS";FP$ 

1340  OPEN  FF7  FOR  t^jT  AS  #1 
1350  N-l :R(0 , D-IPTNRT-DT 
1360  INPUT  “J ,R(N,2) 

1370  IF  R(N,2)  —  1  THEN  CLOSE  #l:GOTO  1490 
1380  R(N,l)~R(N*l,l;-iDT 
1290  N-N+l:GOrO  1360 


1400  PRINT  "EN1TR  THE  PRESSURE  DATA  IN  PSI ,  TERMINATE  INPUT  WITH  A  NEGATIVE" 
1410  PRINT  "VALUE  FOR  THE  PRESSURE," 

1420  PRINT  : PRINT 

1430  N=1 :R(0 , 1)=TSTART-DT 

1440  PRINT  "ENTER  THE  PRESSURE  FOR  POINT  NUMBER  ";N 
1450  INPUT  L(N,2) 

1460  IF  R(N, 2)<0  THEN  1490 
1470  R(N,1)~R(N-1,1)+DT 
1480  E-W+l : GOTO  1440 
1490  N-N-I 

1500  I.PRIN1  CHS$(12) 

1510  LPRTNT  "  HIGH  PRESSURE  COMBUSTION  RESEARCH  TEAM" 
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1520  LPRINT  "  BALLISTIC  RESEARCH  LABORATORY" 

1530  LPRINT  "  INTERIOR  BALLISTICS  DIVISION  -  ABC  BRANCH" 

1540  LPRINT  "  MINI  CLOSED  BOMB  INFORMATION  AND  RESULTS" 

1550  LPRINT : LPRINT 

1560  LPRINT  TAP '40) ; "IGNITER" ; TAB (60) ; "PROPELLANT" 

1570  LPRINT  "If  ITUS(FT-LB/LB) :  "  ; TAB (37)  ; 

1580  LPRINT  USI  '  »###### . #####" ; FI/12 ; : LPRINT  TAB(60);"NOT  COMPUTED" 
1590  LPRINT  "FLAME  TEMP . (K) : " ; TAB (37) ; 

1600  LPRINT  USING  "###### . #####" ; TI ; : LPRINT  TAB(57); 

: LPRINT  USING  "######.#####"; TP 
1610  LPRINT  "MOLECULAR  WEIGHT: " ;TAB(37) ; 

1620  LPRINT  USING  "######.#####"; MI; : LPRINT  TAB(57); 

: LPRINT  USING  "######.#####" ;MP 
1630  LPRINT  "CO-VOLUM£(INCHA3/LB):";TAB(37); 

1640  LPRlNT  USING  "###### . #####" ; NI ; : LPRINT  TAB(57); 

: LPRINT  USING  "######.#####" ;NP 
1650  LPRINT  "WEIGHT: ";TAB(37); 

1660  LPRINT  USING  "###### . #####" ; Cl ; : LPRINT  TAB(57); 

: LPRINT  USING  "######.#####" ; CP 

1670  LPRINT  "DENSITY(LB/INCHA3) :";TAB(37) ;"N0T  COMPUTED" ; TAB (57) ; 

1680  LPRINT  USING  "######.#####" ;RHO 
1690  LPRINT: LPRINT 

1700  LPRINT  TAB ( 20 ); "GRAIN  INFORMATION ( INCHES ) " 

1710  LPRINT  "DIAMETER: "; TAB (18) ; 

1720  LPRINT  USING  "######.#####"; D 
1730  LPRINT  "LENGTH:"; TAB (18); 

1740  LPRINT  USING  "###### .  #####"  ;  L 
1750  LPRINT  -'NUMBER  OF  PERFS"  ;  TAB (18)  ; 

1760  LPRINT  USING  "######.#####" ;NNP 
1770  LPRINT  "AVG.  PERF  DIAM: " ;TAB(18) ; 

1780  LPRINT  USING  "###### .  #####"  ;  PD 
1790  LPRINT  "OUTER  WEB: " ; TAB (18) ; 

1800  LPRINT  USING  "###### . #####" ; WO 
1810  LPRINT  "MIDDLE  WEB: " ;TAB(18) ; 

1820  LPRINT  USING  "######.#####" ;WM 
1830  LPRINT  "INNER  WEB: " ; IAB(18) ; 

1840  LPRINT  USING  "######.»####" ;WI 
1850  LPRINT  "FORM  FUNCTION: " ;TAB(18) ; 

1860  LPRINT  USING  "######. #*###” ;FFO 
1870  LPRINT : LPRINT 

1880  LPRINT  TAB (20) ; "BOMB  AND  PRESSURE  INFORMATION" 

1890  LPRINT  "BOMB  VOLUME(CC) : " ,TAB(22) ; 

1900  LPRINT  USING  "######.#####" ;VB 
1910  LPRINT  "OBSERVED  PMAX(PSI) : " ;TAB(22) ; 

1920  LPRINT  USING  »###### . #####" ; PMAX 
1930  LPRINT  "THEO..  PMAX(PSI)  : "  ;TAB(22)  ; 
t 940  LPRINT  USING  "###### . #####" ; PTHEO 
1950  REM  **  SMOOTHING  PRESSURE  DATA  SO  THE  PRESSURE 
IS  ALWAYS  INCREASING  ** 

1960  REh  *■*** ******************************************************** 
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1970  FOR  1=1  TO  N-l 

1980  IF  R(I , 2)<-R(I+l ,2)  THEN  2010 

1990  IF  I-N-l  THEN  R(I+1 , 2)=R(I , 2) :GOTO  2010 

2000  IF  R(I,2)<»  R(I+2 ,2)  THEN  R(I+1 , 2)-(R(I , 2)+R(I+2 , 2) )/2  ELSE 
R(I,2)-(R(I-l,2)+R(I+l,2))/2 
2010  NEXT  I 

2020  REM  **  COMPUTATION  OF  PROPELLANT  BURNT  AT  EACH  TIME  STEP  ** 

2030  REM  ******************************************************************* 
2040  FOR  1=1  TO  N 

2050  R(I,4)=(R(I,2)*CONl-CON2)/(a,J3+R(I,2)*CON4) 

2060  NEXT 

2070  CLS: PRINT  : PRINT  : PRINT  : PRINT 

2080  PRINT  "DO  YOU  WANT  A  FIVE  POINT  SMOOTHING  BRIDGE  FOR  dP/dt  AND  dM/dt?" 
2090  INPUT  "ANSWER  Y  OR  N";T$ 

2100  IF  T$0"Y"  THEN  2240 

2110  REM  **  COMPUTING  DP/DT  AND  DC/DT  USING  A  FIVE  POINT  BRIDGE  ** 

2120  REM  ********************************************************* 

2130  DEF  FNF(1)=( -2*R(I -2 ,Q1) -R(I-1 ,Q1)+R(I+1 ,Ql)+2*R(I+2 ,Ql) )/(10*DT) 

2140  FOR  Ql»2  TO  4  STEP  2 
2150  FOR  1=3  TO  N-2 
2160  R(I,Q1+1)=FNF(I) 

2170  NEXT  I 
2180  NEXT  Q1 

2190  R(2,3)=(R(2, 2) -R(l, 2) ) /DT:R(N-1 , 3)=(R(N-1 , 2) -R(N-2 . 2) )/DT 

2200  R(N, 3)=(R(N, 2) -R(N-1 , 2) )/DT:R(2 , 5)='R (2 ,4) -R (1 ,4) )/DT 

2210  R(N-l,5)=(R(N-l,4)-R(N-2,4))/DT:R(N,5i=-(’K.(N,4)  -R(N-1 ,4) )/DT 

2220  INPUT  "ANOTHER  SMOOTHING  PASS?";Z$:IF  Z$="Y"  THEN  2140 

2230  GOTO  2280 

2240  FOR  1=2  TO  N 

2250  R(I,3W(R(I,2)-R(I-1,2))/DT 

2260  R(I,5)=(R(I,4)-R(I-1,4) )/DT 

2270  NEXT  I 

2280  REM  **  COMPUTATION  OF  SURFACE  AREA  BASED  ON  FORM  FUNCTION  ** 

2290  REM  *********************************** k***^: **********************>:*** 

2300  IF  FFO-1  THEN  AREA=ATN ( 1 ) *D*D : FOR  1=1  TO  N:R(I, 9)=AREA:NEXT  I 

2310  IF  FFO-1  THEN  NOG=CP/(AREA*L*RHO) ;GCTO  2360 

2320  IF  FFO-2  THEN  GOSUB  2860 

2330  IF  FFO-3  THEN  GOSUB  2960 

2340  T.F  FF0=4  THEN  GOSUB  3070 

2350  IF  FFO-5  THEN  GOSUB  3260 

2360  REM  **  FINDING  AVERAGE  TIME  AND  PRESSURE  VALUES  ** 

2370  REM  ************************ ***i *************************-x***********x 
2380  FOR  1=2  TO  N 

2390  R(I,?)=((R(I,1)+R(I-1.1))/2)*1000 
2400  R(I , 8)=(R(I , 2)+R(I- 1 , 2) )/2 
2410  NEXT 

2420  REM  COMPUTATION  OF  BURNING  RATE  ** 

2430  FOR  1=2  TO  N 

2440  R( I , 6 ) =R ( I , 5 ) /R ( I , 9 ) /RHO/NOG 

2450  NEXT  I 


57 


l-t  trm  1 1 u  ■  if  ■  i'*  t*va' "X  mn rxr 


i  ■u«u»tivu« 


2460  CLS: PRINT  : PRINT  : PRINT  : PRINT 

2470  PRINT  11  ENTER  A  NAME  FOR  THE  OUTPUT  FILE." 

2480  INPUT  "DO  NOT  INCLUDE  A  DRIVE,  DATA  WILL  BE  SAVED  ON  THE  A  DRIVE" ;N$ 
2490  OPEN  "A: "+N$  FOR  OUTPUT  AS  #2 

2500  CLS : LPRINT  "NUMBER  OF  GRAINS  " ,NOG : LPRINT  CHR$(12) 

2510  PRINT  "TIME" ."PRESSURE" , "DP/DT %" , "R" , "SURF.  A" 

2520  PRINT  #2 , "TIME" , "PRESSURE" , "DP/DT" , "  %" , "R" , "SURF.  A" 

2530  PRINT  "MSEC" , "PSI" , "PSI/SEC" , "BURNT" , "IN/SEC" , "FRAC. " 

2540  PRINT  #2, "MSEC" , "PSI" , "PSI/SEC" , "BURNT" , "IN/SEC" , "FRAC. " 

2550  FOR  1=2  TO  N 

2560  PER=INT( (R(I ,4) /CP)*10000)/100 

2570  PRINT  R(I , 7) , INT(R(I ,8)),R(I,3), PER,R(I , 6) ,R(I , 9)/R(l , 9) 

2580  PRINT  #2 ,R(I , 7) , INT(R(I ,8))1R(I,3), PER,R(I , 6) ,R(I , 9)/R(l , 9) 

2590  NEXT  I 

2600  PRINT  #2,CHR$(12) :CLS 

2610  PRINT  "WOULD  YOU  LIKE  A  TABLE  OF  PRESSURE  VS  RATE  STARTING  AT  THE  "; 
2620  PRINT  INT(R(2,8)/1000)*1000+1000; "  PSI  INCREMENTED  BY  500  PSI?" 

2630  INPUT  "ENTER  Y  OR  N";A$ 

2640  IF  A$="N"  THEN  2850 
2650  PRINT  : PRINT  : PRINT 
2660  LPRINT  "PRESSURE" , "RATE" 

2670  PRINT  #2, "PRESSURE", "RATE" 

2680  LPRINT  "PSI" , "IN/SEC" 

2690  PRINT  #2 , "PSI" , "IN/SEC" 

2700  CON9=INT(R(2 , 8)/1000)*1000+1000 
2710  CON10=INT(R(N, 8)/1000)*1000 
2720  C0N11=R(N, 8) -CONIO 

2730  IF  CONll>=500  THEN  CONI 2=CON10+500  ELSE  CON12-CON10 

2740  FOR  I=CON9  TO  C0N12  STEP  500 

2750  FOR  J=2  TO  N 

2760  IF  I<=R(J,8)  THEN  2780 

2770  NEXT  J 

2780  IF  J=N  THEN  2810 

2790  RATE=( (I-R(J-1 , 8) )/(R(J , 8) -R(J-1 , 8) ) )*(R(J , 6) -R(J-1 , 6) )+R(J-l, 6) 

2800  GOTO  2820 
2810  RATE  =R(N, 6) 

2820  LPRINT  I, RATE 
2830  PRINT  #2, 1, RATE 
2840  NEXT  I 
2850  CLOSE  #2: RETURN 

2860  REM  **  SUBROUTINE  TO  COMPU”F  SURFACE  AREA  FOR  SINGLE  ** 

2870  REM  **  PERF  INHIBITED  ON  ALL  SIDES  EXCEPT  THE  PERF  ** 

2880  AREA=3.14159*PD*L 

2890  VOL=3. 14159* (D*D-PD*PD)*L/4 

2900  NOG=CP/VOL/RHO 

2910  FOR  1=1  TO  N 

2920  TCZ=R(I ,4) /NOG 

2930  R(I , 9)=L*3 . 14159*(PD*PD+4*TCZ/(L*RHO*3 . 14159) ) A . 5 
2940  NEXT  I 
2950  RETURN 
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2960  REM  SUBROUTINE  FOR  SINGLE  PERF 
2970  ?I=4*ATN(1) 

2980  V0Lr=PI*L/4*(D*D-PD*PD) 

2990  NOG=C P /V OL/RHO 

3000  FOR  1=1  TO  N 

3010  TCZ=R(I,4;/NOG 

3020  BB=( PD*PD - D*D - 2*L*PD - 2*L*D) 

3030  TEMP=( -BB- SQR(BB*BB-4*(4*D+4*PD)*(2*TCZ/(RHO*PI) ) ) )/(8*(D+PD) ) 
3040  R(I,9)=(PI/2)*((D-2*TEMP)A2-(PD+2*TEMP)A2)+PI*(L-2*TEMP)*(PD+D) 
3050  NEXT  I 
3060  RETURN 

3070  REM  SUBROUTINE  FOR  N  PERF  WITH  NOl 
3080  PI=4*ATN(1) 

3090  V0L=PI*L/4*(D*D-NNP*FD*PD) 

3100  N0C-CP/RH0/70L 
3110  FOR  1=1  TO  N 
3120  TCZ=R(I ,4) /NOG 

3130  ALPHA- (L*NNP -L-2*D- 2*NNP*PD ) / ( 2 - 2*NNP) 

3140  BETA— (D*D-NNP*PD*PD+2*NNP*L*PD+2*L*D)/(4-4*NNP) 

3150  GAMMA— TC2/(RH0*PI*(2-2*NN?) ) 

3160  TE11P1-0 

3170  TEMP-TEMPI- ( TEMP ]. A 3+ALPHA*TEMPl A 2-rBETA*TEMPl+GAMMA)  / 

( 3*TEMP1 a2+2*ALPHA*TEMP1+BETA) 

3180  IF  ABS(T£MP-TEMP1)<. 000001  THEN  3200 
3190  TEMP1-TEMP: PRINT  TEMP: GOTO  3170 
3200  PRINT  "THE  VALUE  OF  R  IS  " ; TEMP 
3210  GOTO  3230 

3220  REM  FINDING  SURFACE  AREA 

3230  R(I,9)=(PI/2)*((D-2*TEMP)A2-NNP*(PD+2*TEMP)A2)+PI*(L-2*TEMP) 

* ( NNP* ( PD+2*TEMP ) + ( D - 2*TEMP ) ) 

3240  NEXT  I 
3250  RETURN 

3260  REM  SUBROUTINE  FOR  SURFACE  AREA  OF  INHIBITED  SINGLE  PERF 
3270  REM  WHERE  INHIBITION  IS  ON  OUTER  LATERAL  SURFACE  ONLY 
3280  PI=4*ATN(1) 

3290  V0L=PI*I./4*(D*D-PD*PD) 

3300  NOG-CP/VOL/RHO 
3310  FOR  1=1  TO  N 
3320  TCZ=R(I,4)/N0G 
3330  ALPHA- (-L/2+PD) 

3340  BETA=( -L*PD/2 -D*D/4+PD*PD/4) 

3350  GAMMA-TCZ/ (RHO*PI*2 ) 

3360  TEMP1-0 

3370  TEMP-TEMPI - (TEMPI A 3+ALPHA*TEMPl A 2+BETA*TEMPlfGAMMA) 

/ ( 3*TEMP1 A  2+2*ALPHA*TEMPl+BETA) 

3380  IF  ABS(TEMP-TEMP1)<. 000001  THEN  3400 
3390  TEMP1-TEMP: PRINT  TEMP.-GOTO  3370 
3400  PRINT  "THE  VALUE  OF  R  IS  " ; TEMP 
3410  GOTO  3430 

3420  REM  FINDING  SURFACE  AREA 
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3430  R ( I , 9 ) -PI* ( PD+2*TEMP)* ( L- 2*TEMP) + ( FI/2 ) * ( D*D - ( PD+2*TEMP ) * ( PD+2*TEMP ) ) 
3440  NEXT  I 
3450  RETURN 
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VARIABLE  LIST 


ALPHA 

AREA 

A$ 

BB 

BETA 

Cl 

CP 

CONI 

CON2 

CON3 

C0N4 

D 

D 

FFO 

FI 

FP 

FP$ 

GAMMA 

KK 

L 

MI 

MP 

N 


Temporary  Storage 

Initial  Surface  Area  of  One  Grain 

Indicator  For  a  TAble  of  Burning 
Rates  by  Pressure 

Temporary  Storage 

Temporary  Storage 

Weight  of  Igniter 

Weight  of  Propellant 

Stores  a  Constant 

Stores  a  Constant 

Stores  a  Constant 

Stores  a  Constant 

Grain  Diameter 

Time  Increment 

Form  Function  Indicator 

Theoretical  Impetus  of 
Igniter 

Experimental  Impetus  of 
Propellant 

Input  File  Name 
Temporary  Storage 

Loop  Index 

Grain  Length 

Molecular  Weight  of  Igniter 
Molecular  Weight  of  Propellant 
Number  of  Data  Points 
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VARIABLE  LIST  (Con't) 


NI 

NNP 

NOG 

NP 

N$ 

0 

PD 

PI 

PIG 

PMAX 

PTHEO 

RHO 


TCZ 

TEMr 

TEMPI 

TI 


Co-volume  of  Igniter 
Number  of  Perfs 
Number  of  Grains 
Co -volume  of  Propellant 
Name  of  Output  File 
Selected  Option 
Average  Perf  Diameter 
Value  of  pi 

Pressure  Due  to  Igniter 
Observed  Maximum  Pressure 
Theoretical  Maximum  Pressure 


Propellant  Density 


Array  Used  to  Store  Values 
R( - , 1)  --  Time 
R( - , 2)  --  Pressure 
R(- , 3)  --  dP/dt 

R(-,4)  --  Total  Mass  of  Propellant 
Burnt 

R(-,5)  --  dC/dt,  Mass  Generation 
Rate 

R(-  ,6)  --  Burn  Fate 
R(-.7)  --  Average  Time 
R(-,8)  --  Average  Pressure 
R(-,9)  --  Surface  Area 


Total  Mass  Burnt  From  a  Single 
Grain 


Temporary  Storage 
Temporary  Storage 
Flame  Temperature  of  Igniter 


b2 
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VARIABLE  LIST  (Con't) 


TOS 

Weighted  Average  of  Flame 
Temperatures 

TP 

Flame  Temperature  of  Propellant 

TS 

System  Temperature 

TSTART 

Starting  Time 

T$ 

Indicator  of  Using  Five  Point 
Bridge 

VB 

Empty  Bomb  Volume 

VOL 

Initial  Volume  of  a  Single  Grain 

WI 

Inner  Web 

WM 

Middle  Web 

WO 

Outer  Web 

-$ 

Mode  of  Data  Input 
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APPENDIX  C 

SAMPLE  OUTPUT  FROM  MINICB 


TABLE  C-l.  Example  of  output  from  MINICB  using  the 
experimental  impetus  option. 


IMPETUS  CALCULATION 


IGNITER  IMPETUS: 
IGNITER  CO- VOLUME: 
IGNITER  WEIGHT: 


BOMB  VOLUME: 
MAXIMUM  PRESSURE: 


390962  FT-LB/LB 
26.15  CU  IN/LB 
.002205  LB 

.05292  LB 
26.15  CU  IN/LB 

100  CC 
55479  PSI 


PROPELLANT  WEIGHT: 
PROPELLANT  CO -VOLUME: 


IMPETUS-  1168.419  J/G 
IMPETUS-  390896  FT-LB/LB 


TABLE  C-2.  Example  of  disk  file  output  from  MINICB 
with  burning  rate  option  chosen. 


TIME 

PRESSURE 

DP/DT 

% 

R 

SURF. 

MSEC 

PSI 

PSI/SEC 

BURNT 

IN/SEC 

FRAC. 

.1 

108 

65000 

.54 

4.707604 

1 

.3 

122 

75000 

.62 

5.431525 

1 

.5 

138 

80000 

.69 

5.793249 

1 

.7000001 

155 

89999.99 

.78 

6.516937 

1 

.9000001 

174 

99999.99 

.88 

7.240458 

1 

1.1 

195 

110000 

.98 

7.9638 

1 

1.3 

218 

120000 

1.1 

8.686944 

1 

1.5 

243 

130000 

1.22 

9.40985 

1 

1.7 

270 

145000 

1.36 

10.49439 

1 

1.9 

301 

160000 

1.52 

11.57852 

1 

2.1 

334 

175000 

1.69 

12.66222 

1 

2.3 

371 

lq5C00 

1.88 

14.10712 

1 

2.5 

412 

210000 

2.08 

15.1897 

l 

2.7 

455 

225000 

2.3 

16.27169 

1 

2.9 

503 

25C000 

2.54 

18.07603 

1 

3.100001 

555 

270000 

2.8 

19.51781 

1 

3.300001 

611 

295000 

3.08 

21.31992 

1 

3.500001 

675 

325000 

3.39 

23.48195 

1 

3.700001 

740 

345000 

3.73 

24.91989 

1 

3.900001 

813 

380000 

4.09 

27.43965 

1 
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TABLE  C.-2.  Example  of  disk  file  output  from  MINICB 
with  burning  rate  option  chosen.  (Con't) 


4.100001 

392 

410000 

4.49 

29.59606 

4 . 30000i 

977 

445000 

4.92 

32,11097 

4.500002 

1069 

473000 

5.37 

34.26246 

4.700001 

1169 

520000 

5.87 

37.49264 

4.900002 

1277 

560000 

6.41 

40.35829 

5.ioocor* 

1393 

605000 

6.99 

43.58008 

l . 300002 

1519 

650000 

7.61 

46.79668 

5 .  ‘'00002 

1654 

700000 

8.29 

50.36783 

5.700002 

1799 

755000 

9.01 

54.29203 

5 . 900002 

1956 

810000 

9.79 

58.20872 

6 . 100003 

2124 

874999  9 

10.62 

62.83522 

6 . 300002 

2305 

934999.9 

11.32 

67.09289 

6 . 500002- 

2500 

1010000 

12.08 

72.4153 

6.70C003 

,'.’.708 

1075000 

13.31 

77.0082 

6.900003 

2932 

1164000 

14.62 

83.01926 

7.100C03 

3171 

1235000 

15.8 

38.?.Qei6 

7.300004 

3423 

1330000 

17  06 

94 , 93803 

7.50U003 

3703 

1420000 

18. U 

101.2970 

7.700004 

399? 

1520000 

10.fi  6 

108.2989 

7.900004 

■oil 

1620000 

21.39 

115.2722 

8.100003 

4546 

373*000 

23.04 

123.281b 

8.300GC3 

5005 

1885000 

24.79 

131.6106 

8.500002 

5388 

1-375000 

76.66 

139.9005 

8.700G02 

5797 

21150C0 

28.65 

149.5617 

8.900001 

6234 

2250000 

30.77 

158.8186 

9 . 100001 

6699 

2400000 

33.02 

169.0787 

9.3C0001 

7195 

2560000 

35.42 

179.9785 

9.500001 

7724 

2730000 

37.98 

191,50bl 

9.700001 

3>87 

2905000 

40.65 

203.3073 

9.899999 

8887 

3095000 

43.57 

216. 0C6 

10.1 

9526 

3290000 

46.62 

229.0725 

10.3 

10206 

3510000 

49,67 

243.703 

10.5 

10930 

3730000 

?3.21 

258.2045 

10.7 

11700 

3570000 

56.97 

273.9456 

10.9 

12519 

4225000 

60.84 

250.5591 

11.1 

13391 

M90090 

64.95 

308.0206 

11.3 

14319 

4780000 

69  3 

326.3C53 

11.5 

15305 

5085000 

73.91001 

342. 726 

11.7 

16354 

5405000 

73.79 

365.9094 

11.9 

17470 

5750000 

83.95 

38V. 4975 

12.1 

18657 

6120000 

89 . 43 

410.4453 

12.3 

19919 

6505000 

95.21 

434.0344 

12.5 

21110 

5400000 

100 

358.5719 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

f, 

1 

1 

\ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

] 

1 

1 

1 
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TABLE  C-2.  Example  of  disk  file  output  from  MlNsCS 
with  burning  rate  option  chosen  (Oon't) 


PRESSURE  RATE 

PSI  IN/SEC 


1000 

32.63715 

1500 

46.30971 

2000 

59.41682 

2500 

72.4153 

3000 

84.51807 

3500 

96.6403 

4000 

108.3656 

4500 

119.7843 

5000 

131.483 

5500 

142.5343 

6000 

153.8562 

6500 

164.6878 

7000 

175.6933 

7500 

186.626 

8000 

197.2873 

8500 

207.826 

9000 

218.3577 

9500 

228.5429 

10000 

239.2708 

1G500 

249 . 3917 

110GC 

259.6352 

Ii500 

269.  ?r7 

12000 

2G0.0274 

12590 

290.1638 

3  3000 

300.1809 

13200 

310.1595 

14000 

320.0166 

14500 

329.8686 

15000 

3.39.7118 

15500 

349 ,4o83 

16000 

359.0886 

16500 

3G8.7252 

1V000 

378.4017 

17500 

388.0775 

18000 

397.7438 

18500 

407.41 

19000 

416.854 

19500 

426.1963 

20000 

358.5719 

2H500 

358.5719 

21000 

358.5719 
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TABLE  C-3.;  Example  of  hardcopy  output  from  MINICB 
using  burning  rate  option. 


HIGH  PRESSURE  COMBUSTION  RESEARCH  TEAM 
BALLISTIC  RESEARCH  LABORATORY 
INTERIOR  BALLISTICS  DIVISION  -  ABC  BRANCH 
MINI  CLOSED  BOMB  INFORMATION  AND  RESULTS 


IMPETUS (FT -LB/LB) : 
FLAME  TEMP.(K): 
MOLECULAR  WEIGHT: 

CO -VOLUME ( INCH*  3 /LB) : 
WEIGHT: 

DENSITY(LB/INCHA3) 


IGNITER 
224583.30000 
2000.00000 
24,77100 
30.00000 
0.00000 
NOT  COMPUTED 


PROPELLANT 
NOT  COMPUTED 
3092.00000 
24.77100 
26.62000 
0.01650 
0.05600 


DIAMETER: 

LENGTH: 

NUMBER  OF  PERFS 
AVG.  PERF  DIAM. 
OUTER  WEB: 
MIDDLE  WEB: 
INNER  WEB: 

FORM  FUNCTION: 


GRAIN  INFORMATION (INCHES) 
0.50000 
1.50000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
1.00000 


BOMB  VOLUME (CC): 
OBSERVED  PMAX(PSI): 
THEO.  PMAX(PSI) : 
NUMBER  OF  GRAINS 


BOMB  AND  PRESSURE  INFORMATION 
59.00000 
21650.00000 
21743.79000 

1.000367 


PRESSURE 

PSI 

1000 

1500 

2000 

2500 

3000 

350U 

4000 

4500 

5000 


RATE 

IN/SEC 

32.63715 

46.30971 

59.41682 

72.4153 

84.51807 

96.6403 

108.3656 

119.7843 

131.483 
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5500 

6000 

6500 

7000 

7500 

8000 

8500 

0000 

9500 

10000 

10500 

11000 

11500 

12000 

12500 

13000 

13500 

14000 

14500 

15000 

15500 

16000 

16500 

17000 

17500 

18000 

18500 

19000 

19500 

20000 

20500 

21C00 


TABLE  C-3.  Example  of  hardcopy  output  from  MINICB 
using  burning  rate  option.  (Con't) 


142.5343 

153.8562 

164. 68',  8 

175.6933 

186.626 

197.2873 

207.826 

218.3577 

228.5429 

239.2708 

249.5917 

259.6355 

269.857 

280.0274 

290.1638 

300.1809 

310.1595 

320.0166 

329.8686 

339.7118 

349.4683 

359.0886 

368.7252 

378.4017 

388.0775 

397.7438 

407.41 

416.354 

426.1961 

358.5719 

358.5719 

358.5719 
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APPENDIX  D 


TRADITIONAL  MEASUREMENT  OF  EXPERIMENTAL  IMPETUS 


In  this  approach,  heat  loss  and  effects  of  the  igniter  are  ignored. 
Problems  with  the  igniter  can  be  eliminated  by  having  the  igniter  the  same  as 
the  propellant.  The  following  symbols  will  be  used: 


P 

V 

F 

n 

L 

c 


Pressure 

Chamber  Free  Volume 
Impetus 
Co -volume 
Loading  Density 
Charge  Weight 


Start  with  the  Noble-Able  equation  of  state, 
P(V  -  nc)  -  Fc 


(D-l) 


Divide  by  c, 


P(V/c  -  n)  -  F 


(D-2 


Next  divide  by  P  to  obtain 


V/c.  -  n  =  (1/P)  F 


(D-3) 


Since  loading  density  is  the  charge  weight  divided  by  volume,  L  *=  c/V. 
Substituting  into  equation  (D-3)  gives 
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1/L  -  Fa/P)  +  n 


(D-4) 


Now  by  performing  a  series  of  closed  bomb  firings  with  different  loading 
densities  and  recording  the  maximum  observed  pressure  a  plot  of  1/L  versus  1/P 
can  be  made.  The  plot  should  be  linear  and  the  slope  will  be  the  impetus  and 
the  intercept  the  co -volume. 
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APPENDIX  E 


EQUATIONS  FOR  THE  HAND  CALCULATOR  VERSION  OF  MINICB 


As  mentioned  in  the  main  body  of  the  report  an  alternate  approach  to 
computing  the  burning  rate  is  to  determine  the  distance  burnt,  x,  which 
corresponds  to  each  pressure  and  then  compute  dx/dt  directly  instead  of  using 
equation  (17),  r  -  (dCp/dt)/(p*A) .  For  the  hand  calculator  version  of  MINICB 
this  approach  was  utilized  since  it  requires  fewer  programming  steps.  The 
analysis  is  identical  to  the  point  of  finding  Cp  ..  equation  (14) ,  for  each 
value  of  pressure.  Next,  it  is  necessary  to  determine  values  for  x  for  each 
value  of  Cp.  In  the  analysis  presented  earlier,  values  of  x  were  found  by 
deriving  a  relation  between  x  and  Cp  for  the  specific  propellant  under 
investigation  based  upon  initial  grain  geometry  and  the  relation 

Cp  =  density  *  Volume  Used  (E-l) 

Unfortunately,  solving  the  expression  resulting  from  equation  E-l  for  x  was, 
for  most  grain  geometries,  quite  complicated,  often  requiring  a  numerical 
scheme.  The  number  of  programming  steps  required  for  such  an  operation  on  a 
hand  calculator  was  felt  to  be  excessive  and  would  result  in  a  separate 
program  for  each  grain  geometry.  Thus,  it  was  decided  to  make  a  further 
simplification  in  the  analysis  for  the  hand  calculator  version  so  as  to  allow 
for  the  widest  possible  range  of  form  functions.  Starting  with  equation  (E- 
1) ,  divide  by  the  original  mass  of  propellant  giving  equation  (E-2). 


Cp/C  -  density  *  Volume  Used/C  (E-2) 

Now  Cp/C  is  che  mass  fraction  burnt  and  will  be  denoted  by  M.  Also,  C,  the 
total  mass  of  propellant,  is  the  density  multiplied  by  the  total  original 
volume.  Thus,  equation  (E-2)  can  be  rewritten  to  produce  uquaticn  (E-3)  .. 
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M  “ 


density  *  Volume  Used 


(E-3) 


density  *  Total  Original  Volume 

Canceling  density  and  writing  the  volumes  in  terms  of  volume  per  grain 
equation  (E-4)  can  be  obtained,. 

Number  of  Grains  *  Volume/grain  used 

M  -  -  (E-4) 

Number  of  Grains  *  Original  Volume/grain 


herefore . 

M  *  Original  Volume/grain  -  Volume/grain  used  (E-5) 

Now  the  original  volume/grain  is  a  constant  based  upon  the  original  grain 
geometry  and  volume/grain  used  is  determined  by  the  original  grain  geometry 
and  the  distance  burnt,  x.  Hence,  equation  (E-5)  could  be  used  to  define  a 
relation  between  M,  the  mass  fraction  burnt,  and  the  depth  burnt,  x.  Using 
this  expression  to  actually  calculate  x  would  require  as  many  programming 
steps  as  before.  However,  assuming  increasing  values  of  x,  a  reasonable 
assumption,  a  series  of  values  for  the  mass  fraction  burnt  can  be  computed. 
These  values  may  then  be  fitted  empirically  to  an  equation  such  as 


Biawiii  Mff’aflgwc:aBBMWTtgi<wBWgi»BawBfB^grw»»w 


TrtPfrv 


Routines  for  such  a  curve  fitting  are  available  for  many  calculators .  The 
empirical  equation  can  then  be  used  to  determine  the  values  of  x  which 

correspond  to  each  value  of  Cp  in  the  program.  This  approach  will  reduce  the 

computation  of  x  to  an  exponentiation  and  multiplication.  Routines  for  such  a 
curve  fitting  are  available  for  many  calculators.  Results  for  the  seven  perf 

grain  used  in  validation  5,  Table  5  are  shown  in  Table  E-l.  Again  the 

analys{s  is  only  valid  up  to  the  point  of  slivering  which  for  a  seven  perf 
grain  occurs:  at  about  80%  all  burnt. 


TABLE  E-l.  Fitting  to  an  Equation 

x  ->  a*M3  b*M2  +  c*M  +  d 
for  a  Se^en  Perf  Grain 


Web  -  0.03988  inches 


Distance  Burnt  Mass  Fraction  Burnt 


Inches 

.0025 

.09195951 

.005 

.1883845 

.0075 

.2891035 

.01 

.3939452 

.0125 

.502738 

.015 

.6153106 

.0175 

.7314915 

.01994 

.8481996 
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Values  for  a,  b,  c,  and  d 


a  -  0.001758388 
b  -  -0.006253527 
c  -  0.0275523105 
d  -  0.00002226^ 

Goodness  of  Fit:  P.^  -  0.999999885 


Using  the  technique  above  for  various  other  form  functions  has  resulted 
ivi  similar  valued  for  R  .  Tavs,  it  appears  that  using  a  cubic  equation  to 
relate  depth  burnt  and  mass  fraction  burnt  is  a  reasonable  assumption  and 
should  provide  an  acceptable  alternative  for  determining  the  depth  burnt,  x, 
in  the  hand  calculator  version  of  MINICB. 


APPENDIX  F 

HAND  CALCULATOR  PROGRAM 
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Listings  for  the  various  programs  and  subprograms  iniich  make  up  the  hand 
calculator  version  of  MINICB  are  given  below  for  a  Hewlett-Packard  HP-41C 
calculator.  The  remaining  section  of  the  Appendix  presents  comparisons 
between  runs  of  the  hand  calculator  version  and  the  BASIC  version  of  MINICB 
and  a  comparison  with  CBRED2. 


Hewlett-Packard  and  HP-41C  are  tradenames  of  the  Hewlett-Packard 
Corporation.  Use  of  these  names  does  not  constitute  indorsement  by  the  US 
Army . 


PROGRAM  LISTING  --  HP-41.C 


Procedure  Name 

Function 

TINYCB 

Calls  Various  sub¬ 
programs 

ENDATA 

Enter  Data 

CONST 

Computes  All  Constants 

ENPRES 

Enter  Pressure 

PVSRle 

Computes  Rates 

PVSR2 

Interpolation  of  Rates 

MINICB 

(Hand  calculator  version) 


0 i ♦18L 

"TIHVCB" 

82 

XEQ 

•EHBRT8- 

8; 

m 

“COHST” 

84 

XEQ 

"EHPRES* 

85 

XE? 

-pvjPie- 

86 

9 

87 

EMI? 

85 


ENDATA 


8!*LSL  ’EHDflTA* 

92  129 
9:  STO  91 

04  (4y 
85  S’ 3  92 
96  168 
8?  ST0  94 
98  183 
09  STO  93 

19  -IGNITER  BflTfl.,.- 

11  flVIEK  ' 

12  -IMPETUS’- 

13  PROMPT 

14  STO  81 

15  -fl.  TEMP, 7 * 

19  PROMPT 

17  SVO  02 

18  -CCVOLUPE?* 

19  PROMPT 
29  STO  93 

21  'MOL.  HT.?* 

22  PROMPT 

27  STO  84 

2*  ’HEIGHT?" 

25  DRGHFT 

28  STO  05 

27  'PROPELLANT  Dfi’P" 

28  PVIEH  - 

2?  ’DENSITY?’ 

38  °RQMPT 


31  SVC-  80 
?2  ’FL.-THHP.-” 

33  °ROHPT 

34  STO  97^ 

35  ’COyOLUME?’ 

36  PROMPT 

37  STO  88 

33  -MOL.  HT.'1’ 
3S  PROMPT 

48  STO  89 

41  "WEIGHT?" 

42  PROMPT 

43  STO  18 

4i  ’  GROIN  'j-vfi. 

45  flVIEfi 

46  ’BIPHETER?’ 

47  PROHPT 
43  STO  11 

49  •LENGTH?" 

59  PROMPT 

51  STO  12 

52  " GENERAL  3RTP 

53  flVIEU 

54  -OSS.  PMhX?* 

55  PROMPT 

56  ST0  13 

57  ’BOMB  ¥&.*• 

58  PPQflPT 

59  STO  14 
68  EHt 


CONST 


9l*LSL  "C0Hjt* 

46  - 

82  12 

4?  370 

1 “l 

87  ST*  Si 

58  PCL 

1« 

8-  73372 

51  is.; 

3786 

85  PCL  87 

52  / 

36  * 

53  Rt  l 

10 

87  PCL  18 

54  RCL 

96 

0S  •' 

55  • 

8?  ;CL  89 

56  - 

I'3  ' 

57  RCL 

95 

cr1  a: 

ca  n 

wC-  K'.l 

i  «*  r  •*  ’ 

5C  * 

17  * 

88  - 

<4  - 

61  370 

77 

.b 

1:  K _  1- 

52  3C. 

95 

*6  ‘  6 . 7i78'; 

63  RCL 

8: 

*  7 

64  * 

IS  ?C'_  1C 

65  RCL 

70 

l?  PCL  8? 

66  * 

28  4 

67  RCL 

92 

o;  - 

63  / 

22  RCL  85 

69  3TQ 

34 

23  RCL  93 

78  RCL 

83 

24  t 

71  RCL 

86 

25  - 

72  !A' 

26  •• 

77  - 

27  270  79 

74  3T0 

36 

2?  ?CL  :8 

75  RCL 

34 

2?  RCL  37 

76  RCL 

33 

39  - 

77  •' 

31  PCL  95 

73  3T9 

7"y 

.» 

i  ■*  r  i  .■} 

*C 

7?  PCI 

[J 

7^  > 

38  :(t> 

7d  - 

8*  ?’ 

?c:  19 

32, 

3t  »C’-  95 

C7  i 
w-. 

77  i 

84  ■ 

»,  w 

35  270 

7C 

J  “  *■ 1  -1  !•  i 

36  pr>_ 

12 

40  3CL  13 

q  ;  * 

41  » 

36  'CL 

0b 

42  5CL  -30 

$2  at 

w  j 

?8 

i '  77 

**■**■  v  >  -J  Ji- 

0!  P." 

.  •  ■  -t  w 

*5 

«=*  777"  * 

*♦.' 

92  - 

lli,  f 

97  j’G 

7’ 

47  ;r .  jc 

■3*i  «:n 

87 


ENPRES 


91+lSL  '£::P5E:'’ 

?£  ‘fit? 

87  ;'0  39 
0J  *B£:-Th  Tl^E-" 
0:  3C0fJC’ 

00  3T7  16 

r  1800 

86  67-  It 
9°  "MO.  OF  PT“ 
18  PF0FP7 

11  '-TO  2-3 

12  6T0  28 

1-j  «88i 

14  37*  28 

-  C  \ 

14  ST*  28 

t  -1  n  -  . 

l>  if 

!§  ;7 3  21 
1?  2CL  37 

29  *ps,  1  >* 

21  SKL  X 

K  Py,]HpT 

23+13L  30 

24  RCL  21 

25  X 

26  'VI  3 

27  ‘PRE5.  HC.  * 
2"  iF.CL  X 

2V  PP0P°T 

30  Fin  3 

31  370  IBB  38 


iO'J 

»7C 

EHL 


39 

21 

30 


88 


r.w 


u  unu  jsmavwumr  Kff, 


PVSRle 

01 *181  *PYCRle* 

82  8 

83  STO  48 
04  4 

85  STO  41 
06  100 
07  STO  99 
88  *1  CONST?* 

09  PROMPT 
18  STO  42 

11  *X  CONST?* 

12  PROMPT 

13  STO  43 

14  *Xt2  CONST?- 

15  PROMPT 

16  STO  44 

17  *Xt3  CONST?* 

18  PROMPT 

19  STO  45 

28  *  .PR.  .RATE. * 

21  PROMPT 

22  RCL  20 

23  STO  21 
24*L8L  48 

25  RCL  IHB  98 

26  RCL  33 

27  * 

28  RCL  34 

29  - 

30  RCL  INS  98 

31  RCL  36 

32  * 

33  RCL  35 

34  < 

35  / 

36  RCL  18 

37  / 

38  STO  INS  91 

39  RCL  INS  91 
48  3 

41  YfX 

42  RCL  45 

43  * 

44  RCL  INC  91 

45  Xt2 

46  RCL  44 

47  * 

48  f 

49  RCL  INC  91 

50  RCL  43 

51  * 


52  * 

53  RCL  42 

54  + 

55  STO  INC  92 

56  RCL  21 
5?  1.9 

58  X>Y? 

59  GTO  41 

58  RCL  INS  92 

61  I 

62  ST-  92 

53  RCL  INS  92 

64  CHS 

65  RCL  2 

66  + 

67  RCL  16 

68  / 

63  STO  IHB  94 

70  RCL  INC  98 

71  1 

72  ST-  90 

73  RCL  INC  90 

74  RCL  Z 

75  + 

76  2 

77  / 

78  STO  INC  93 

79  SLA 

89  RCL  INC  93 

81  FIX  INC  49 

82  (RCL  X 

83  *F  * 

34  RCL  IKS  94 
85  FIX  INC  41 
36  ARCL  X 
87  PROMPT 


88 

1 

89 

ST+ 

38 

90 

ST+ 

92 

91+LBL 

41 

92 

1 

93 

ST+ 

98 

94 

ST+ 

91 

95 

v  V 

ST+ 

92 

96 

ST+ 

93 

97 

ST+ 

94 

98 

ISO 

21 

99 

GTO 

40 

188 

FIX 

9 

181 

END 

89 


INTERPOLATION  PROGRAM 
PVSR2 


81 *151  *PVSR2' 

58  * 

00  •  L  « 

82  18! 

5!  S79  64 

183  RCL  56 

8?  STQ  93 

52*LBL  58 

161  “In  4 

84  RCL  INS  93 

53  RCL  23 

162  9PCL  X 

05  5880 

54  1808 

133  5R0flPT 

06  / 

55  / 

18“  ISl'  64 

67  irr 

56  2 

105  rc  58 

88  5680 

57  * 

106  “i:-:  ? 

89  * 

53  $70  22 

137+-.5L  ‘ROTE" 

18  5069 

59*L3L  51 

13?  1 

11  * 

66  183 

109  ?3 

12  TO  53 

6;  STO  ?: 

i  :0  lotf 

13  RCL  23  •' 

62  RCL  22 

i.l  "1  22 

4  ‘  ‘ 

63  INT 

112  IT-7 

>  2“  _ 

64  $7+  93 

1  <7  * 

18  ST+  93 

65  RCL  INS  93 

;  >  S“0  ?- 

I?  RCL  INI  ?7 

66  1808 

ll:  ■'CL  64 

18  5808 

5?  f 

115  ur 

19  /  4"  "  ’ 

26  INT  f 

2!  5899* 

68  RCL  64 

I 1 7  1099 

6?  INT 

11?  * 

76  X<=Y’ 

119  RCi  IHE  33 

22  ♦ 

7f  {■  TA  c-c?  . 

I  11  U.'U  Mi. 

128  - 

23  5880] 

72  ISG  22 

12!  CTO  65 

24  i-  ‘I  . 

'  73  £70  51 

122  RCL  INB  93 

25  S7G  61 

74+L3L  52 

123  CHS 

25  CH5  ' 

;c  pr; 

<  «  t'wu  w 

12“  1 

27  RCL  IHL  93 

76  INT 

125  Ft  93 

28  *  \ 

?**  t 
«  *  * 

126  RCw  IKE  93 

2°  ST0.62 

70  - 

;  0”  p.-  7 

I_]  Viw  L 

36  2580 

70  0Q»  07 

'OC  ^ 

3!  '*.< -V. 

60  X*V? 

*  r  .  •  o  i 

32  070  55 

81  3™t  53 

13?  :  a 

33  RCL  61 

82  XEQ  "PP7E- 

:“!  RCL  6-5 

34  $T0  63 

33  510  54 

*  c2  x 

35  uT0  56 

34*L8L  53 

177  670  6-7 

36+L8L  55 

35  5CL  27 

134  RCL  IhT  94 

37  RCL  61 

86  :5? 

i  . 

35  2560 

87  ♦ 

‘  t  “,T—  '3  4 

39  * 

33  370  3- 

17“  “'Cl  INS  ?4 

40  isl  63 

39  RCL  !l‘D  ?4 

i  7-^  *  ■  -• 

4UL3L  56 

93  ? -0  66 

i jc  sr,_  2 

4-5  RCL  57 

5;  +i  Q‘  "u 

4  JU  * 

43  58 

32  RC.  64 

141  'CL  67 

dJ. 

93  IN’ 

14;  * 

45  1090008 

?i  1000 

1  A;  RCL  JHP  34 

46  ' 

95  1 

■  ill 

A7  RCL  60 

'46  rT>  ? 

!“c  ?“:■  66 

48  i960 

97  CL; 

i  ir  ^ 

4? 

?-*  n".*-  -« 

U7  “j-i:- 
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To  evaluate  the  hand  calculator  version  of  MINICB  three  different  cases 
were  selected.  The  first  two  involved  a  comparison  with  results  from  the 
BASIC  version  of  MINICB.  These  comparison  were  performed  to  investigate  the 
effect  of  expressing  the  depth  burnt  in  terms  of  the  mass  fraction  burnt  via  a 
cubic  function.  Thus,  the  first  comparison  was  with  the  the  cord  burning  in  a 
cigarette  fashion.  This  formunction  results  in  a  perfect  fit  between  x  and 
M,  a=b=d=0  and  c  =  grain  length,  and  hence  there  should  be  almost  perfect 
agreement  between  the  results.  The  only  difference  is  the  numerical  accuracy 
between  the  machines.  Results  for  this  comparison  are  shown  in  Table  F-l  and 
Figure  F-l.  In  the  second  comparison,  a  seven  perf  grain  was  used  to  evaluate 
the  effects  of  the  curve  fitting.  Details  of  the  curve  fitting  for  this  grain 
are  given  in  Appendix  E.  Results  are  found  in  Table  F-2  and  Figure  F-2. 
Finally,  a  comparison  was  made  with  results  of  CBRED2  for  a  19  perf  grain. 

The  results  are  given  in  Table  F-3  and  Figure  F-3.  In  viewing  the  results  it 
is  important  to  keep  in  mind  that  the  hand  calculator  version  was  being 
performed  with  a  very  limited  number  of  pressure- time  points,  11,  24,  and  13 
points  respectively  for  the  comparisons. 


TABLE  F-l.  Comparison  of  Calculator  Version  of 
MINICB  with  BASIC  Version  for  a  Cord 
Burning  in  Cigarette  Fashion 


PRESSURE  RATE  %  Difference 


psi 

in/sec 

MINICB 

MINICB 

(BASIC) 

(CALC) 

2000 

58.90745 

58.9075 

-0.00008 

4000 

107.5714 

107.5713 

0.000092 

6000 

152.6348 

152.6348 

0 

8000 

195.8801 

195.8801 

0 

10000 

237.7167 

237.7167 

0 

12000 

278.4922 

278.4922 

0 

14000 

318.1077 

318.1077 

0 

16000 

357.5178 

357.5178 

0 

18000 

243.6403 

243.1748 

0.191426 

11  data  points 


Paroant  arror 


for  Clf.  grain  —  11  Points 


Figure  F-l.  Comparison  of  Calculator  Version  of 
MINICB  with  BASIC  Version  for  a  Cord 
Burning  in  Cigarette  Fashion 

As  an  absolute  error  of  0.02%  indicates,  the  calculator  version  of  MINICB 
and  MINICB  produce  essentially  the  same  results  for  the  cigarette  burning  form 
function. 
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TABLE  F-2.  Comparison  of  Calculator  Version  of 

MINICB  with  BASIC  Version  for  a  Seven 
Perf  Grain 


PRESSURE  RATE 


%  Difference 


psi  in/sec 


5000 

MINICB 

(BASIC) 

1.967996 

10000 

3.567153 

15000 

5.055268 

20000 

6.464363 

25000 

7.824321 

30000 

9.142696 

35000 

10.4264 

40000 

11.6847 

45000 

12.91465 

50000 

14.12269 

55000 

15.31281 

60000 

16.48633 

65000 

17.63615 

70000 

18.78596 

75000 

19.9091 

80000 

21.01668 

85000 

22.12426 

90000 

23.23184 

95000 

24.29764 

100000 

25.3589 

105000 

26.42017 

110000 

27.48144 

115000 

28.5427 

120000 

29.37056 

MINICB 

(CALC) 

1.9513 

0.857 

3.5548 

0.347 

5.0577 

-0.049 

6.4880 

-0.364 

7.8724 

-0.611 

9.2167 

-0.803 

10.526 

-0.952 

11.811 

-1.068 

13.066 

-1.156 

14.297 

-1.222 

15.510 

-1.274 

16.706 

-1.315 

17.877 

-1.345 

19.047 

-1.371 

20.191 

-1.396 

21.319 

-1.419 

22.448 

-1.440 

23.576 

-1.459 

24.668 

-1.501 

25.756 

-1.542 

26.844 

-1.58 

27.932 

-1.614 

29.021 

-1.647 

29.88 

-1.706 

24  data  points 
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?«ro«nt  trror  of  rt to* 


i 

0.8 
0.6 
0.4 

o.e 

0 

-0.8 
-0.4 
-0.6 
-0.8 

-1.8 
-1.4 
—1.6 
-1.8 

5000  15000  85000  35000  45000  55000  65000  75000  85000  95000  105000  115000 

Praaaora  (pal) 


for  7  Parf  grain  -  84  Data  points 


Figure  F-2.  Comparison  of  Calculator  Version  of 

MTNTCB  with  BASIC  Version  for  a  Seven 
Perf  Grain 


An  absolute  error  of  1.17%  indicates  that  using  a  cubic  function  to  fit  depth 
burnt  and  mass  fraction  burnt  is  a  fairly  reasonable  assumption  for  multi-perf 
grains . 


94 


TABLE  F-3.  Comparison  of  Calculator  Version  of 
MIHICB  with  CBRED2  for  a  Nineteen 
Perf  Grain 


PRESSURE  RATE  %  Difference 


psi 

in/se 

MINICB 

c 

CBRED2 

5000 

1,385 

1.378 

-0.533 

10000 

2.554 

2.543 

-0.432 

15000 

3.723 

3.742 

0.505 

20000 

4.865 

4.875 

0.216 

25000 

6.008 

6.026 

0.304 

30000 

7.147 

7 . 184 

0.57 

25000 

8.283 

8 . 302 

0.232 

40000 

9.411 

9.352 

-0.628 

45000 

10.539 

10.377 

-1.537 

50000 

11.294 

11.45 

.381 

13  data  points  used 


for  19  Porf  groin  -  13  Data  points 


I  Figure  F-3.  Comparison  of  Calculator  Version  of 

I  M.INICB  wj  ch  CBRED?.  for  a  Nineteen 

!  Perf  Grain 

! 

I 

i 

With  an  absolute  error  of  0.63%  the  agreement  between  the  calculator 
version  of  MIN1CB  and  CBRED2  is  excellent  fcr  this  particular  comparison 
funilar  results  have  been  obtained  for  other  form  functions  and  propellants. 

i 

|  Based  upon  the  result  of  the  comparisons  it  is  felt  that  the  calculator 

version  of  KINICB  can  be  used  with  a  high  degree  of  confidence  for  computing 
j  burning  rates  from  a  liraited  number  of  data  points. 
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APPENDIX  G 

COMPARISON  OF  THE  BASIC  VERSION  OF  MINICB 
WITH  CBRED2 
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For  the  sake  of  completeness,  the  results  of  MINICB  were  compared  for  a 
limited  number  of  cases  with  the  results  from  BRL's  standard  burning  rate 
reduction  code,  CBRED2,  an  extension  of  CBP.ED.  Typical  results  are  presented 
below.  A  common  experimental  pressure-time  history  obtained  from  a  closed 
bomb  firing  was  used  as  input  to  both  programs.  Thus,  this  comparison  is 
relative,  comparing  only  the  performance  of  both  programs  to  each  other  and 
not  to  a  burning  rate  known  a  priori.  Details  of  the  propellant  and  igniter 
used  in  the  closed  bomb  firing  are  given  in  Table  G-l.  Results  from  the  code 
are  presented  in  Table  G-2  and  graphically  in  Figure  G-l. 


TABLE  G-l.  Details  of  Propellant  and  Igniter 
used  in  Closed  Bomb  Firing 


Propellant 

Igniter 

Type: 

JA2 

FFFG 

Weight  (lb) 

0.1528 

0.002205 

Density(lb/inA3) : 

0.0578 

Gamma: 

1.2247 

1.25 

Impetus(ft-lb/lb) : 

363241 

97500 

Co-volume(inA3/lb) 

27.42 

30.00 

Flame  Temperature (K) : 

3437 

2000 

Number  of  data  points  101 


TABLE  G-2.  Comparison  of  Burning  Rate  Between 

MINICB  and  CBRED2  for  a  19-perf  Grain 


Pressure 

Rate 

Percent  Difference 

Psi 

Inch/s 

MINICB 

CBRED2 

5000 

1.335 

1.378 

3.185 

10000 

2.489 

2,543 

2.13 

15000 

3.734 

3.742 

0.21 

20000 

4.833 

4.875 

0.87 

25000 

6.072 

6.026 

-0.75 

30000 

7 .211 

7.184 

-0.37 

35000 

8.339 

8.301 

-1.05 

40000 

9.404 

9.352 

-0.55 

45000 

10.352 

10.377 

0.24 

50000 

11.445 

11.45 

0.044 

55000 

12.337 

12.466 

I  .05 

60000 

12.987 

12.929 

-0.45 
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P*re*at  Dlfftrtno* 


Comparison  Between  MINICB  and  CBRED2 


19  P«p t  Qritn 


Figure  G-l.  Percent  Difference  in  Burning  Rate 
Betveen  MINICB  and  C3RFD2 


As  shown  in  Table  G-2  and  Figure  G-l  the  results  for  CBRED2  and  MINICB 
are  in  excellent  agreement  with  an  absolute  average  error  of  about  0.9%. 
Similar  results  were  obtained  for  comparisons  involving  different  propellants 
and  other  form  functions. 
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No.  Of 
Copies 


DISTRIBUTION  LIST 


No.  Of 

Organization  Copies  Organization 


1 2  Commander 

Defense  Technics].  Info  Center 
ATTN:  DTJ.C-L'DA 

Cameron  Station 
Alexandria.-  VA  22304-6145 

1  Commander 

USA  Concepts  Analysis  Agency 
ATTN:  D.  Hardiscn 
8120  Woodmont  Avenue 
Bethesda,  MD  20014 

1  HQDA/DAMA-ZA 

Washington,  DC  20310 

1  HQDA,  DAMA-CSM,  E-  Lippi 
Washington,  DC  20310 

1  HQDA/DAMA-ART-M 

Washington,  DC  203 "SO 

1  HQDA/SARDA 

Washington,  DC  20310 

1  Commander 

US  Army  War  College 
ATTN:  Library-FF2  29 
Carlisle  Barracks,  PA  17013 

1  Director 

US  Army  BMD 

Advanced  Technology  Center 
P.  C.  Box  1500 
Huntsville,  AL  35807 

1  Chairman 

DOD  Explosives  Safety  Board 
Room  856-C 
Hoffman  Bldg.  1 
2461  Eisenhower  Avenue 
Alexandria,  VA  22331 


1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

5  Project  Manager 

Cannon  Artillery  Weapons 
System,  ARDC,  AMCCOM 
ATTN:  AMCPM-CW , 

F .  Menke 
AMCPM-CWW 
AMCPM-CWS 
M.  Fisette 
AMCPM-CWA 
R.  DeKleine 
H.  Hassmann 
Picatinny  Arsenal,  NJ 
07806-5000 

2  Project  Manager 

Munitions  Production  Base 
Modernization  and  Expansion 
ATTN:  AMCPM-PBM,  A.  Siklosi 

AMCPM-P3M-E,  L.  Laibson 
Picatinny  Arsenal,  NJ 
07806-5000 

3  Project  Manager 

Tank  Main  Armament  System 
ATTN:  AMCPM-TMA,  K.  Russell 
AMCPM-TKA- 105 
AMCPM-TMA- 120 
Picatinny  Arsenal ,  NJ 
07806-5000 

1  Commander 

US  Army  Watervliet  Arsenal 
ATTN:  SARWV-RD,  R.  Thierry 
Watervliet,  NY  12189 


1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCPM-GCM-WF 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


DISTRIBUTION  LIST 


No.  Of  No.  Of 

Copies  Organization  Copies  Organization 


22  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-TSS 
SMCAR-TDC 
SMCAR-LC 

LTC  N.  Barron 
SMCAR-LCA 
A.  Beardell 
D.  Downs 
S.  Einstein 
S.  Westley 
S.  Bernstein 
?.  Kemmey 

A.  Bracuti 
J.  Rutkowski 

SMCAR-LCB-I 

D.  Spring 

SMCAR-LCE ,  R. Walker 
SMC.AR-LCM-E 
S.  Xaplowitz 
SMCAR-LCS 
SMCAR-LCU-CT 

E.  Barrieres 
R.  Davitt 

SMCAR-LCU-CV 
C-  Mandala 
SMCAR-LCW-A 
M.  Salsbury 
SMCAR-SCA 
I-.  Stiefel 

B.  Brodman 

Picatinny  Arsenal,  NJ 
07806-5000 

4  Commander 

US  Army  Armament  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 

L.Ambrosini 
AMSKC»IRC,  G.  Cowan 
AMSMC-LEM 
W.  Fortune 
R.  2astrow 

Pock  Island,  JL  51299-7300 


1  Director 

Benet  Weapons  Laboratory 
Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-LCB-TL 

Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Reseat ’h 
and  Development  Command 
ATTN:  AMSAV-E 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120 

1  Commander 

US  Array  TSAP.CCM 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120 

1  Director 

US  Army  Air  Mobility  Research 
And  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

I  Commander 

US  Army  Communications 
Electronics  Command 
ATTN:  AMSEL-ED 
Fort  Monmouth,  NJ  07703 

1  Commander 

ERADCOM  Technical  Library 
ATTN:  STET-L 

Fort  Monmouth,  NJ  07703-5301 
1  Commander 

US  Army  Harry  Diamond  Lab. 
ATTN:  DELHD-TA-L 

2800  Powder  Mill  Road 
Adelphi ,  MD  2078?, 

1  Common der 

US  Array  Missile  Command 
Rsch,  Dev,  &  Engr  Ctr 
ATTN:  AhSMI-IlD 
Redstone  Arsenal,  AL  35893 


DISTRIBUTION  LIST 


No.  Of 

Copies  Organization 

1  Director 

US  Army  Missile  &  Space 
Intelligence  Center 
ATTN:  AIAMS-YDL 
Redstone  Arsenal, 

AL  35898-5500 

1  Commandant 

US  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-TSL 
Warren,  MI  48397-5000 

1  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-CG 
Warren,  MI  48090 

1  Project  Manager 

Improved  TOW  Vehicle 
ATTN:  AMCPM-ITV 
US  Army  Tank  Automotive 
Command 

Warren,  MI  48090 

1  Program  Manager 

Ml  Abrams  Tank  System 
ATTN :  AMCPM-GMC- S A , 

T.  Dear 

Warren,  MI  48090 

1  Project  Manager 

Fighting  Vehicle  Systems 
ATTN:  AMCPM-FVS 
Warren,  MI  48090 

1  President 

US  Tunny  Armor  &  Engineer 
Board 

ATTN:  ATZK-AD-S 
Fort  Knox,  KY  40121 


No.  Of 

’'Copies  Organization 

1  Project  Manager 

M-60  Tank  Development 
ATTN :  AMCPM-M6  0  TD 
Warren,  MI  48090 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range, 
NM  88002 

1  Commander 

US  Army  Training  &  Doctrine 
Command 

ATTN:  ATCD-MA/  MAJ  Williams 
Fort  Monroe,  VA  23651 

1  Commander 

US  Army  Materials 
Technology  Laboratory 
Dyna  East  Corporation 
ATTN:  Christine  P.  Brandt, 
Document  Control 
3132  Market  Street 
Philadelphia,  PA  19104-2855 

1  Commander 

US  Army  Research  Office 
ATTN :  Tech  Library 
P.  O.  Box  12211 
Research  Triangle  Park,  NC 
27709-2211 

1  Commander 

US  Army  Belvoir 
Research  &  Development  Ctr 
ATTN:  STRBE-WC 

Tech  Library  (Vault) 
Bldg  315 

Fort  Belvoir,  VA  22060-5606 

1  Commander 

US  Army  Logistics  Ctr 
Defense  Logistics  Studies 
Fort  Lee,  VA  23801 


103 


DISTRIBUTION  LIST 
No.  Of 

Copied  Organization 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA  62R 
SEA  64 

Washington,  DC  20362-5101 


No.  Of 

Copies  Organization 

1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905 

1  President 

US  Army  Artillery  Board 
Ft.  Sill,  OK  73503 

1  Commandant 

US  Army  Command  and 
General  Staff  College 
Ft  Leavenworth,  KS  66027-5080 

1  Commandant 

US  Army  Special  Warfare 
School 

ATTN:  Rev  &  Tng  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commander 

Radford  Army  Ammo  Plant 
ATTN :  SMCRA-QA/HI  LIB 
Radford,  VA  24141 

1  Commander 

US  Army  Foreign  Science  & 
Technology  Center 
ATTN:  AMXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901 

1  Commandant 

US  Army  Field  Artillery 
Center  &  School 
ATTN:  ATSF-CO-MW,  B.  Willis 
Ft.  Sill,  OK  73503 

1  Commander 

US  Army  Development  and 
Employment  Agency 
ATTN;  MODE-TED-SAB 
Fort  Lewis,  WA  98433 

1  Office  of  Naval  Research 

ATTN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  Street 
Arlington,  VA  22217 


1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-S54-Tech  Lib 
Washington,  DC  20360 

1  Assistant  Secretary  of  the 
Navy  (R,  E,  and  S) 

ATTN :  R.  Reichenbach 
Room  5E787 
Pentagon  Bldg. 

Washington,  DC  20350 

1  Naval  Research  Lab 
Tech  Library 
Washington,  DC  20375 

2  Commander 

US  Naval  Surface  Weapons 
Center 

ATTN:  J.  P.  Consaga 
C.  Gotzmer 

Silver  Spring,  MD  20902-5000 

Naval  Surface  Weapons  Center 
ATTN :  S .  Jacobs/R 1 0 
Code  730 

K.  Kim/Code  R-13 
R.  Bernecker/Code  R-13 
Silver  Spring,  MD  20902-5000 

5  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G33,  J.  L.  East 
W.  Burrell 
J •  Johndrow 
Code  G23,  D.  McClure 
Code  DX-21  Tech  Lib 
Dahlgren,  VA  22448-5000 


104 


DISTRIBUTION  LIST 


No.  Of  No.  Of 

Copies  Organization  Copies  Organization 


2  Commander 
Naval  Underwater 

Systems  Center 
Energy  Conversion  Dept. 

ATTN:  CODE  5B331,  R-  S.  Lazar 
Tech  Lib 

Newport,  RI  02840 

3  Commander 

Naval  Weapons  Center 
ATTN:  Code  388,  R.  L.  Derr 
C.  F.  Price 
T.  Boggs 

Info .  Sci .  Div . 

China  Lake,  CA  93555-6001 

1  Superintendent 

Naval  Postgraduate  School 
Dept,  of  Mechanical 
Engineering 

Code  1424  Library 
Monterey,  CA  93943 

1  Program  Manager 
AFOSR 

Directorate  of  Aerospace 
Sciences 

ATTN:  L.  H.  Caveny 
Bolling  AFB,  DC  20332 

3  Commander 

Naval  Ordnance  Station 
ATTN:  J.  Birkett 
D.  Brooks 
Tech  Library 
Indian  Head,  MD  20640 

1  HQ  AFSC/SDOA 

Andrews  AFB,  MD  20334 

1  AFRPL/DY,  Stop  24 

ATTN:  J.  N.  Levine/DYCR 
Edwards  AFB,  CA  93523-5000 

1  AFRPL/TSTL  (Tech  Library) 

Stop  24 

Edwards  AFB,  CA  93523-5000 


1  AFRPL/MKP3,  Stop  24 
ATTN:  B.  Goshgarian 
Edwards  AF3,  CA  93523-5000 

1  AFFTC 

ATTN:  SSD-Tech  Lib 

Edwards  AFB,  CA  93523 

1  AFATL/DLYV 

ATTN:  George  C.  Crews 
Eglin  AFB,  FL  32542-5000 

1  AFATL/DLJE 

Eglin  AFB,  FL  32542-5000 

1  Air  Force  Armament  Lab. 
AFATL/DLODL 

Eglin  AFB,  FL  32542-5000 
1  AFWL/SUL 

Kirtland  AFB,  NM  87117 

10  Central  Intelligence  Agency 
Office  of  Central  Reference 
Dissemination  Branch 
Room  GE-47  HQS 
Washington,  DC  20502' 

1  Central  Intelligence  Agency 
ATTN:  Joseph  E.  Backofen 
HQ  Room  5F22 
Washington,  DC  20505 

1  General  Applied  Sciences  Lab 
ATTN :  J .  Erdos 
Merrick  &  Stewart  Avenues 
Westbury,  NY  115S0 

1  Aerodyne  Research,  Inc. 
Bedford  Research  Park 
ATTN:  V.  Yousef ian 
Bedford,  MA  01730 

i  Aerojet  Solid  Propulsion  Co. 
ATTN:  P.  Micheli 
Sacramento,  CA  95813 


105 


DISTRIBUTION  LIST 


No.  Of  No.  Of 

Copies  Organization  Copies  Organization 


1  Atlantic  Research  Corporation 
ATTN:  M.  K.  King 
5390  Cheorokoe  Avenue 
Alexandria,  VA  22314 

1  A VCD  Everett  Rsch  Lab 
ATTN:  D.  Stickler 
2385  Revere  Beach  Parkway 
Everett,  MA  02149 

1  Calspan  Corporation 
ATTN:  Tech  Library 
P.  O.  Box  400 
Buffalo,  MY  14225 

1  General  Electric  Company 
Armament  Systems  Dept . 

ATTN:  M.  J.  Bulman j 
Room  1311 
Lakeside  Avenue 
Burlington,  VT  05401 

1  IITRI 

ATTN:  M.  J.  Klein 
10  W.  35th  Street 
Chicago,  IL  60616 

1  Hercules  Powder  Co. 

Allegany  Ballistics 
Laboratory 
ATTN:  R.  3.  Miller 
P.  O.  Box  210 
Cumberland,  MD  21501 

1  Hercules,  Inc 
Bacchus  Works 

ATTN:  K.  P.  McCarty 
P.  0.  Box  98 
Magna,  UT  84044 

2  Director 
Lawrence  Livermore 

National  Laboratory 
ATTN:  M.  S.  L-355, 

A .  Buckingham 
P.  O.  Box  808 
Livermore,  CA  94550 


1  Lawrence  Livermore 

National  Laboratory 
ATTN:  M.  S.  L-355 
M.  Finger 
P.  0.  Box  803 
Livermore,  CA  94550 

1  Olin  Corporation 

Badger  Army  Ammunition  Plant 
ATTN:  R.  J.  Thiede 
Baraboo,  WI  53913 

1  Olin  Corp/Smokeless  Powder- 
Operations 
R&D  Library 
ATTN:  V.  McDonald 
P.O.  Box  222 
St.  Marks,  FL  32355 

1  Paul  Gough  Associates,  Inc. 
ATTN:  P.  S.  Gough 

P.  O.  Box  1614, 

1048  South  St. 

Portsmouth,  NH  03801 

1  Physics  International  Company 
ATTN :  Library 

H.  Wayne  Wampler 
2700  Merced  Street 
San  Leandro,  CA  94577 

2  Rockwell  International 
Rocketdyne  Division 

ATTN:  BAC8  J.  E.  Flanagan 
J .  Gray 

6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

1  Princeton  Combustion  Research 
Lab • ,  Inc . 

ATTN:  M.  Summerfield 

475  US  Highway  One 

Monmouth  Junction,  NJ  08852 


106 


No.  Of 
Copies 


DISTRIBUTION  LIST 


Organization 

1  Science  Applications,  Inc. 
ATTN:  R.  B.  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

3  Thiokol  Corporation 
Huntsville  Division 
ATTN:  D.  Flanigan 
R.  Click 
Tech  Library 
Huntsville,  AL  35807 

1  Scientific  Research 
Assoc. , Inc. 

ATTN:  H.  McDonald 
P.O.  Box  498 
Glastonbury,  CT  06033 

1  Thiokol  Corporation 
Wasatch  Division 
ATTN:  Tech  Library 
n .  O.  Box  524 
Brigham  City,  UT  84302 

2  Thiokol  Corporation 

Elkton  Division 
ATTN :  R.  Biddle 

Tech  Lib. 

P.  O.  Box  241 
Elkton,  MD  21921 

1  Universal  Propulsion  Company 
ATTN:  H.  J.  McSpadden 
Black  Canyon  Stage  1 

Box  1140 

Phoenix,  AZ  85029 

2  United  Technologies 
Chemical  Systems  Division 
ATTN :  R.  Brown 

Tech  Library 
P.  0.  Box  358 
Sunnyvale,  CA  94086 

1  Veritay  Technology,  Inc. 

4845  Millersport  Hwy. 

P.  O.  Box  305 

East  Amherst,  NY  14051-0305 


No.  Of 

Copies  Organization 

1  Battelle  Memorial  Institute 
ATTN:  Tech  Library 
505  King  Avenue 
Columbus,  OH  43201 

1  Brigham  Young  University 

Dept,  of  Chemical  Engineering 
ATTN:  M.  Beckstead 
Provo,  UT  84601 

1  California  Institute  of  Tech 
204  Karman  Lab 
Main  Stop  301-46 
ATTN:  F.  E.  C.  Culick 
1201  E.  California  Street 
Pasadena,  CA  91109 

1  California  Institute  of  Tech 
Jet  Propulsion  Laboratory 
ATTN:  L.  D.  Strand 
4800  Oak  Grove  Drive 
Pasadena,  CA  91103 

1  Professor  Herman  Krier 
Dept  of  Mech/Indust  Engr 
University  of  Illinois 
144  MEB;  1206  N.  Green  St. 
Urbana,  IL  61801 

1  University  of  Minnesota 
Dept,  of  Mechanical 
Engineering 
ATTN:  E.  Fletcher 
Minneapolis,  MN  55455 

1  University  of  Massachusetts 
Dept,  of  Mechanical 
Engineering 
ATTN :  K.  Jakus 
Amherst,  MA  01002 

1  Case  Western  Reserve 
University 

Division  of  Aerospace 
Sciences 
ATTN :  J .  Tien 
Cleveland,  OH  44135 


No.  Of 
Copies 


DISTRIBUTION  LIST 


No.  Of 

Organization  Copies  Organization 


3  Georgia  Institute  of  Tech 
School  of  Aerospace  Eng. 

ATTN:  B.  T.  Zinn 
E  Price 
W.  C.  Strahle 
Atlanta,  GA  30332 

1  Institute  of  Gas  Technology 
ATTN:  D.  Gidaspow 
3424  S.  State  Street 
Chicago,  IL  60616 

1  Johns  Hopkins  University 
Applied  Physics  Laboratory 
Chemical  Propulsion 
Information  Agency 
ATTN:  T.  Christian 
Johns  Hopkins  Road 
Laurel,  MD  20707 

1  Massachusetts  Institute  of 
Technology 

Dept  of  Mechanical  Engineering 
ATTN :  T .  Tonng 
77  Massachetts  Avenue 
Cambridge,  MA  02139 

1  Pennsylvania  State  University 
Dept.  Of  Mechanical 
Engineering 
ATTN:  K.  Kuo 
University  Park,  PA  16802 

1  University  of  Michigan 
Gas  Dynamics  Lab 
Aerospace  Engr  Bldg 
ATTN:  Dr.  G.  M.  Faet.h 
Ann  Harbor,  MI  48109-2140 

1  Purdue  University 

School  of  Mechanical 
Engineering 
ATTN:  J.  R.  Osborn 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 


1  SRI  International 

Propulsion  Sciences  Division 
ATTN :  Tech  Library 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Rensselaer  Polytechnic  Inst. 
Department  of  Mathematics 
Troy,  NY  12181 

1  Director 

Los  Alamos  National  Lab 
ATTN:  M.  Division,  B.  Craig 
T-3  MS  B216 
Los  Alamos,  NM  87545 

1  Stevens  Institute  of 
Technology 
Davidson  Laboratory 
ATTN:  R.  McAlevy,  III 
Castle  Point  Station 
Hoboken,  NJ  07030 

1  Putgers  University 

Dept-  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  S.  Temkin 
University  Heights  Campus 
New  Brunswick,  NJ  08903 

1  University  of  Southern 

California 

Mechanical  Engineering  Dept. 
ATTN:  0HE200,  M.  Gerstein 

Los  Angeles,  CA  90007 

2  University  of  Utah 

Dept,  of  Chemical  Engineering 
ATTN:  A.  Baer 

G.  Flandro 

Salt  Lake  City,  UT  84112 

1  Washington  State  University 
Dept,  of  Mechanical 
Engineering 
ATTN:  C.  T.  Crowe 
Pullman,  WA  99163 


108 


No.  Of 
Copies 


DISTRIBUTION  LIST 


Organization 


No.  Of 

Copies  Organization 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
Cdr,  USATECOM 

ATTN:  AMSTE-TO-F 

Cdr,  USACSTA 

ATTN:  S.  Walton 
G.  Rice 
D.  Lacey 
C.  Herud 

Dir,  HEL 

ATTN:  J.  Weisz 
Cdr,  CRDC,  AMCCOM 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ _ Date  of  Report 

2.  Date  Report  Received _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 
other  area  of  interest  for  which  the  report  will  be  used.) 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

°LD  Organization 

ADDRESS 


Address 


City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


